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This paper investigates the preliminary ship design 
process. The thesis begins by looking at methodologies 
for ship design and concludes that a different outlook on 
the overall process could provide an improvement. In par- 
ticular, possible modifications were identified in the 
o^vner's requirements as traditionally given to the designer. 
After the new design procedure was developed, effort 
shifted to implementing the procedure for a containership 
design. 

A preliminary ship design model is developed. This 
model v/as then used in a test design problem in which an 
owner desires a single ship to add to an existing trade 
route. The design model was used to identify acceptable 
ship alternatives. 

The determination of attractive designs was accom- 
plished by a version of an optimization method known as 
a Sequential Unconstrained Minimization Technique. This 
method employs a barrier penalty function to handle the 
problem constraints. The search method was found to be 
effective in identifying design alternatives. In parti- 
cular, it behaved well with the discontinuities imposed 
by the discrete engine selections involved with gas tur- 
bine power plants. 
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INTRODUCTION 



Through the years, ocean transport systems have provi- 
ded a common forum in v/hich the businessmen and engineers 
could participate as a team. The complexity of the present 
systems as v/ell as the sophistication of their design has 
not evolved without numerous growing pains. Over the years, 
society has developed the framework in v/hich each party 
operates. The boundaries of responsibility for the two are 
not always well defined nor are their liabilities. This 
vagueness extends to their interaction with the rest of soc- 
iety. The results of past efforts to formalize their rules 
with respect to the ship design process can be found in 
rules of the American Board of Shipping and of many insur- 
ance companies. These rules guide the engineer/designer 
in his task so that the final ship design will meet certain 
minimum standards. These standards have been developed over 
the years and are based on experience. 

This paper does not propose a change in the traditional 
roles of the businessman and the engineer as related to 
ocean transport systems, nor does it develop any new innova- 
tions into either of their professional areas. It does, 
hov/ever, attempt to define explicitly the formal relation- 
ships betv/een the parties involved in the design process and 

to indicate hov/ these relationships can be used in actual 
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ship design. The approach taken is not traditional. The 
result is not a nev/ relationship, but rather a more defined 
and process oriented approach where there can be a better 
agreement as to each parties role with the result of improv- 
ing their joint effort. 

In an attempt to demonstrate this philosophy, the the- 
sis presents a sample problem. With this problem, a method 
of solution is then pursued. The important aspects of this 
problem are identified in the context developed earlier in 
the paper which defines the roles of the participants. The 
problem is approached in an objective manner and the outputs 
are identified which would prove useful to a decision maker. 
The paper presents the results with caution. First, because 
the solution addresses only the problem outlined in the 
scenario. Secondly, the design model has not yet been veri- 
fied by present design practice. Finally, the optimization 
method may not in itself identify the overall best alter- 
native to a real life problem. For this reason, no attempt 
is made to identify such a result. 

One of the easiest ways to present a problem is to re- 
late a scenario. From this, the assumptions and simplifi- 
cations made in the analysis become more palatable and the 
solution is not regarded as simply a mass of asstimptions 
presented for the convenience of the solution technique. The 
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scenario also identifies a clearer division of roles of the 
businessman (referred to as the customer) and the engineer 
(referred to as the designer) in the scenario. By restrict- 
ing the problem as stated in the scenario, the real life 
complications and ambiguities are minimized. However, these 
will not be ignored in the development of the design process. 
The information flows of the design methodology guided the 
development of the design program. This program utilized 
the capabilities of the computer and the solution technique. 
Before the design methodology is outlined or the scenario 
developed, it is important to identify the approach of the 
paper. 

This paper uses a systems analysis format for presenta- 
tion, By this it is meant that the first issues covered 
involve the problem definition and objectives of the study. 
Following the problem definition, the paper proceeds to 
identify criteria and measures of effectiveness for the 
problem solution. This is then followed by an outline of 
the different types of alternatives available. From this 
position, a decision model is developed. Tools known as 
optimization techniques are applied to this model. The re- 
sults and the analysis are accomplished by the application 
of a Sequential Unconstrained Minimization Technique re- 
ferred to as SUIvlT. The results are then studied and pre- 
sented so that a decision can be made. This decision may 
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not cover the problems originally raised in the scenario, 
but should at least pave the way for the next design itera- 
tion. As such, the paper demonstrates the computer's capa- 
bilities, This, in turn, will demonstrate the possible role 
of a computer in a computer aided design. 

Definitions are imperative for a proper understanding 
of any study, A list of model variables have been collected 
with their definition in Table found in Appendix A. Their 
use in context should not prove difficult for the reader. 
There is one definitbn, however, that desemres explanation 
before we begin. This relates to the use of the terms 
'closed form’ or 'explicit' as qualifiers of expressions in 
the mathematical model. These will be used interchangeably 
in the paper. By an explicit expression, it is meant that 
the relationships betv/een the variables involved can be 
found in terms of an algebraic expression. In particular, 
we will be interested in the relationships involving the de- 
cision variables. 

The complexity of real life design problems made it 
difficult to identify and assess all the implications of 
the proposed design methodology. The results of the design 
model v/ere useful in identifying the usefulness in a con- 
tainer ship design. This paper should serve as a basis 
from which to reanalyze the existing methods used in the 
design of ship systems. 
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CHAPTER I. 

OCEAN TRANSPORTATION SYSTEMS 



Much has been written about transportation systems. 

The important features of any transportation system are its 
collection, warehousing, transport and distribution subsys- 
tems as well as the associated management and human subsys- 
tems. In today's economy, some transportation systems are 
so large that neither a single individual nor a single cor- 
poration can control and direct the action of all the ele- 
ments. This is especially true of the ocean transportation 
systems. V/e will limit our discussion to the actual trans- 
port by the ocean vehicle. The character of this ocean 
transport subsystem which distinguishes it from other modes 
of travel is the need to have a transporting vessel capable 
of translating an ocean. In this transit, a ship will ex- 
perience a hostile environment. The advantage of ocean 
transport as opposed to other transport methods is the abil- 
ity of an ocean transport system to carry larger loads, 
either in bulk or v/eight. 

1 . 1 Elements 

From the following figure, the ocean transportation sys- 
tem can be divided into three basic parts, “The first is the 
mode of ocean transport by ship. At each end of the ocean 

transport is a port and terminal system. At these terminal 
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Figure 1 OCEAN TRANSPORTATION SYSTEI-1 
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facilities, the cargo is handled so that another mode of 
travel can be utilized. This constitutes the second ele- 
ment. In the last element, we will lump all of the remain- 
ing parts of a general transportation system. These will 
include the various distribution and intermediate transpor- 
tation systems,^ 

1 . 2 Subsystem Interactions 

We want to look now at the interactions between the 
three parts. As described, the elements of the transpor- 
tation system function in series, thus the capacity and the 
related efficiency of the total system can be affected by a 
single element. Any system design will attempt to size each 
of the elements so that there will be no disparity in capa- 
city between single units because of the resulting ineffi- 
ciencies. An exception would be the planning of excess 
capacity to be utilized during a future period of growth. 
Conceptually this model should prove useful and as more 
ports and shipping routes are opened, the same principle 
can be used to describe the larger total system. This sim- 
ple model ignores the management problem of capacity utili- 
zation in large networks, but the interactions should be more 
easily understood. 

• The major purpose of the ship is to act as a mode of 
containment and transportation between ports. The terminal 

also serves as a location to inventory cargo so tfiat loading 
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and unloading of the ship can be achieved without delay. 

Many methods have been proposed for reducing the cargo han- 
dling at the terminals and thus increasing their capacity 
and efficiency. One such attempt is the containerization 
of cargo. Instead of handling the different cargoes indivi- 
dually, containers are used, thus simplifying the loading 
problems and reducing breakage and loss of cargo. Over the 
last few years this method has become even more popular, 
especially with customers that transport small, high value 
items, such as electronic equipment. For large bulk cargoes 
ship barge systems have been developed. 

The interaction between the ship and terminal has many 
different facets. They can be related to a time delay in 
port. V/hen a ship arrives, it av/aits space next to the ter- 
minal, Once along side, the crane will continuously handle 
containers during the loading operation. The duration of 
this evolution depends directly upon the reach, carrying 
capacity and speed of the crane. After the incoming cargo 
is off loaded and the outgoing cargo is on loaded, the ship 
may have to v/ait for favorable navigation conditions such as 
slack v/ater at high tide so that the ship can maneuver out 
of port safely. The interaction is completed in the same 
way that it started, with the ship at sea proceeding between 
ports at its operating speed. 

The interactions betv/een the port terminal and the rest 



17 



of the transportation system can be defined along similar 
lines. It becomes evident that there is a large management 
function at the terminals. The skill developed in these 
functions and the capacities of each of the elements deter- 
mines the volume of cargo flov/ in the total system. To say 
that the transportation system is capacity limited by a sin- 
gle element should require that the element is operating 
with maximum efficiency, Otherv/ise, capacity can be in- 
creased without additional capital investment. The complex 
problem of capcity determination will play an important role 
in the total transportation system design. It should be ob- 
vious that when the total system is operating at capacity, 
subsystems with excess capacity are not desirable unless 
there are future plans for expansion. This points out the 
requirement that the design cover the whole of the transpor- 
tation system, Many papers have addressed themselves to 

the overall capcity balancing problem. The reader is re- 

. 2 2 
f erred to Erichsen and Hancock", 

Independent of the overall transportation system design, 
there v/ill be a stage in the total design when the capabili- 
ties and costs are desired of a ship that can interact with 
a given terminal configuration. This may be for an addition 
of shipping capacity on an already existing route or for an 
entirely nevi transportation system. We are ' interested in 
this interaction between the customer and the designer. 
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CHAPTER II. 



SHIP DESIGN PROCESS 



2, 1 General 

The ship design process is that set of actions involved 
from the initial desire to build a ship iintil the final deli- 
very of the product. In some circumstances, one or both of 
these events may be difficult to determine. The process 
involves two groups. One is the owner or user group and the 
other is the designer or builder group. The traditional pro- 
cess (see Figure 2) involves the determination of the owner's 
expectations and his statement as to the capabilities desired 
in the ship design. Many times this v/ould be the first for- 
mal step in the process. The need for the ship may result 
from intuitive feel5.ngs based upon years of experience as a 
ship operator, V/hatever the source, a gap is identified be- 
tween supply and demand for services. The owner's require- 
ments have traditionally been quite explicit — the speed, 
endurance, and payload of the design being specified. This 
statement provides the means by which the design passed to 
the engineer from the owner. As depicted in the block dia- 
gram, the designer then proceeds to implement a design pro- 
cedure that results in a feasible design which meets the 
owner's requirements. At that point, the owner would decide 
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SHIP DESIGN PROCESS 




Figure 2 
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if the ship design shoiald continue to a building phase. 

The steps listed in Figure 2 are representative of the 
process. They may either be found explicitly or implicitly 
in the design effort. There are two factors which disting- 
uish the process. In general, it is a single pass system. 
Inside each phase there may be iteration, but there is no 
provision of changing the inputs to any phase. The second 
characteristic is the parallel design feature of the 
other elements of the system. The decision to build does 
not necessarily provide for a system wide evaluation after 
the owner's requirements are fixed. Before criticizing this 
form, let us investigate the nature of the system, 

2.2 Normative Design Process 

From the ship design process, the decisions that must be 
made are identified. They are in turn associated with a 
level of decision making. The following conceptual model 
exemplifies the hierarchy of decision making that makes up 
the ship design process (See Figure 3 )* According to 
M. D. Mesarovic, et. al., there is no single best model for 
describing the multilevel, hierarchial system. However, 
the essential characteristics are a "...vertical arrangement 
of subsystems, and [a^ dependence of the higher level sub- 

Zj, 

systems upon the actual performance of the lov/er levels.” 

In the Figure, there is no significance to the number 
of levels drawn. They are used only to show possible 
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INFORT.IATION FLOV/S 
IN A NORIvlATIVE DESIGN PROCESS 



DECISION MAKER 



OWNER 



OWNER 



DESIGNER 



DESIGNER 




CRITERIA 



QUALITATIVE 

CRITERIA 



QUALITATIVE 

CRITERIA 



quantitative 

CRITERIA 



Figure 3 
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relationships. The input to the process is the result of an 
analysis which identifies a need. The output of the process 
is a solution (design) which hopefully meets this need to 
the satisfaction of the owner. At the higher levels, the 
owner plays a major role. At the lower levels, the designer 
accomplishes the design procedure. From each level, require- 
ments and guidelines are passed down to the next lower level, 
and from each lower level the current results and tinresolved 
decisions are passed up for action. The levels of decision 
making have been identified as different stages at which 
tradeoffs between criteria are accomplished. Associated 
with each level is a set of criteria which are related to 
the decisions being made. The decisions concerning certain 
criteria can be delegated to lower stages in the design. In 
this model of the design process, the owner has responsibil- 
ity over all decision levels. Those involving quantitative 
critieria such as cost have been identified with the lower 
stages in the process. This assumes that the owner’s pre- 
ferences have been identified. 

In the current design practice, the designer receives 
the owner's requirements and stops as soon as a good solu- 
tion is reached. This means that the last tv/o stages of 
the design process are in equilibrium with each other, but 
not with the rest of the system. That is, the output of the 
designer is a feasible solution which satisfies the original 
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owner's requirements. The traditional process terminates 
here. If, instead of the last stages being conducted inde- 
pendent of the rest of the system, continuous feedback and 
adjustment were possible, a more desired result may be ob- 
tained, This would require more coordination of effort 
between the owner and the designer throughout the design. 

By setting up a formal system to handle the various 
levels of decision making, the owner will help control the 
whole design process by giving better direction to lower 
level management. This would be most important in the de- 
sign of new ship types where the final desired solution is 
not obvious to the owner. 

This suggests that one of the first steps in the design 
process is the determination of a set of criteria and order- 
ing them in a way which reflects their importance to the 
owner. In the normative process, some of these criteria 
are assigned to lower levels of decision making. Usually 
the criteria can be categorized as either quantifiable or 
non-quantif iable. The quantifiable criteria can be easily 
delegated to lov/er levels of decision making due to the 
straightf orvra.rd procedures involved in the tradeoffs. This 
does not prevent some or all of the nonquant if iable criteria 
from also being assigned. The owner can never comfortably 
remove himself from the tradeoffs required by the nonquanti- 
fiable criteria. He may, however, quantify his preferences 
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so that lower levels can perform tradeoffs similar to those 
accomplished for quantifiable criteria in the decision pro- 
cess, An example of such a delegation is in the economic 
criteria which identifies the owner’s preferences for 
changes in the timing of cash flows. By identifying a dis- 
coxmt factor consistant with the owner's time value of 
money, the designer can proceed with the economic evaluation 
as if the criteria v/ere entirely quantifiable. This relieves 
the owner of the need to make tradeoffs of a large number 
of cases but still allows him to influence the process. An 
example of this influence would be a change in the discount 
rate after several iterations. This method of manipulating 
the economic criteria would reflect uncertainty of the owner 
in his actual preferences. 

For an individual as the owner, the decision process 
presents few conflicts of interest. As the number and 
diversity of those in the decision group grows, internal 
conflicts arise when preferences are expressed. This be- 
comes a significant problem when the owner is as large as 
the government. The socio-political environment may then 
prevent an explicit definition of the preferences due to 
the nature of the system being designed. The problems of 
decision making involved in this environment are felt most 
strongly by those who act as representatives for the govern- 
ment. The decisions based on preferences will produce focal 
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points for conflicting interest groups. This tends to re- 
duce the area in which explicit decisions are made to those 
in which mutual agreement on the method of quantifying the 
criteria can be achieved. This often results in the quali- 
tative criteria being either ignored or being considered in- 
formally. 

This multilevel, hierarchial system is compatible with 
the ship design process. It has the advantages of isolating 
the decision maker and setting responsibilities. This then 
leads to reliable subsystems. In particular, the multilevel 
system offers an improvement over other systems in the uti- 
lization of resources when solving large scale complex prob- 
lems, It also provides more flexibility and in general can 
be expected to produce better system output. This does not 
suggest that the system is ideal. Some of the disadvantages 
are its complex operation or behavior. As such, it is diffi- 
cult to analyze or comprehend. It also is difficult for one 
to control or influence its progress. However, the multi- 
level, hierarchial system has proven that it’s general quali- 
ties are better than other systems designed to handle large 
programs . 

2,3 Proposed Design Process 

Any proposed system must provide a method for the flow 
of information from higher authority, performance of its 
design task, and the reporting of results using the developed 
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criteria to the next higher level of decision making, 
upv/ard flow of information v/ill be used to determine if there 
is need for further iteration to achieve a better overall 
design. This is done by investigating tradeoffs betv/een 
various criteria and their associated costs. The informa- 
tion that would be necessary to accomplish this task v/ould 
include first, the values of the decision variable and an 
estimate of the confidence limits for these values. Secon- 
dly, the information should include a sensitivity analysis 
that v/ould show the decision maker the resulting effect on 
the criteria measures for changes in the values of selected 
variables. The proposed design process (Figure attempts 
to incorporate these requirements. 

The best way to understand the major features of the 
proposed change is to contrast this design process with the 
traditional process covered earlier. First, the proposed 
design process, being iterative, takes more time and effort, 
V/here the designers objectives are fairly clear in the 
traditional process, the guidelines nov/ provide for much 
more leev/ay in the actual design. This in turn requires 
additional effort. The process requires increased interac- 
tion between the ov/ner and designers so that the designer 
can determine the preferences of the owner. The ov/ner *s 
identification of the criteria and measures of effectiveness 
should be given substantial attention so that this 
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MODEL OF THE PROPOSED 
SHIP DESIGN PROCESS 




Figure 4 
28 





interaction is effective. 

In both systems the owner has the responsibility of 
making his desires known and making the final decision to 
build the ship. In the traditional process, the owner would 
use a previous ship design to communicate his desires. When 
this luxury no longer exists, it becomes necessary to for- 
malize the method by which such preferences were made. This 
is especially true for large systems when there is no one 
single decision maker who can perform the vital integration 
function. 

One way to rationalize the prevalence of the tradi- 
tional system is to study the time history of ship design. 
Over the years the design process identified successful ship 
types. This meant that the individual design efforts did 
not require formal iteration to achieve good designs. By 
perpetuation and extrapolation of the ship features that 
had proven themselves in actual service, sufficient success 
in the ship deisng effort could be achieved. This identi- 
fies a situation in which the traditional approach may be 
preferred. V/hen the changes in the need as percieved by 
the ov/ner are fairly constant over time, newer ship designs 
can be based on the more successful existing ships. A 
change in the traditional method would thus require that 
new requirements be placed on the ocean transportation 
system. These would be similar to the needs generated which 
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caused the huilding of liquid natural gas (LNG) ships v/here 
no similar ship type previously existed. 

In order that the additional benefits identified by the 
normative model may be incorporated in the ship design pro- 
cess, a new method was proposed. The major differences were 
first, the iterative nature of the process; secondly, the 
change in the form of the ovmer's requirements entering the 
design procedure. This means that the designer will no 
longer identify merely a feasible solution. In turn, this 
requires a change in the perspective of the designer and 
the ov/ner. No change was made in the actual design proce- 
dures. The next chapter identifes a sample program which 
will be used to demonstrate the proposed design process. 
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CHAPTER III. 



PROBLEM STATEMENT 



3. 1 Scenario 

A large U, S. owner, an operator of containerships, 
desires to expand operations between San Diego, California, 
and Yokohama, Japan, The company is interested in purchas- 
ing a single ship to increase the shipping capacity of this 
route. The company presently has no available ships to 
transport the increased demand for cotton and leather goods 
in Japan, If present rates prevail, it is assumed that each 
ton of cargo transported will generate $100 in revenues for 
the company. At the present time, there are no arrangements 
anticipated for the backhaul of cargo. The company usually 
does not consider the backhaul when making their capital in- 
vestment decisions because the quantity of cargo handled on 
the return voyage is so small that almost any reasonably 
sized ship would be capable of providing the capacity. Each 
alternative would then generate the same backhaul revenues, 
and thus would not affect the decision to be made. 

The trade route imposes some restrictions on the alter- 
natives. First, the terminal facilities are such that a 
single loading sequence will not be able to provide more 
than 1200 containers without excessive delays. Secondly, 
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the pierside cranes maximum reach restricts the ships beam 
to 110 feet. Finally, the harbor channel restricts the 
ship's draft to be less than 37*. The ports have an average 
charge associated with the pilotage and tug fees which is 
incurred upon entering or departing a port. Past experience 
estimates these to average $1,000 per port call, V/hile in 
port, there are additional wharfage and port fees which are 
determined by the number of days in port. The port delay 
times have historically averaged two days per port call, 
V/hile in port, there are additional v;harfage and port fees 
which are determined by the number of days in port. The 
port delay times have historically averaged tv;o days per 
port call. 

The container size in use is a 20' X 8' X 8' container 
which meets the ,'lmerican Standards Association requirements. 
Each full container has been estimated to weigh l6,7 tons. 
This restricts the stacking of containers, one on top of 
another v/ithout an intervening, supporting deck. 

In order that all costs and prices be consistant, the 
base year of 19^9 was selected. The prevailing economic 
conditions in the U,S. in 19^9 were such that the operator 
could be assured of obtaining a 55 per cent building sub- 
sidy from the Maritime Administration upon request,-^ In 
addition, the company was anticipating financing three quar- 
ters of the unsubsidized costs by incurring a 25 year loan 
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at ten per cent v/hich would be insured by the government. 
Other alternative routes have been generating a yield on 
the companies original investment after tax of about tv/enty 
per cent. The company feels that any new route should 
generate a similar yield. In addition to the subsidy, the 
company will be able to take a seven per cent investment 
tax credit. "Or comparisons between alternatives, this will 
be translated into a reduction of initial out-of-pocket 
costs. 

3. 2 Desie:n Problem 

The above scenario presents several questions to be 
ansv;ered. Here we will attempt to identify the characteris- 
tics of the design v/hich would be obtained during a preli- 
minary design stage. This effort will require tliat design 
constraints be translated so that non-feasible alternatives 
can be eliminated from consideration. After a choice is 
made of the ship's characteristics, its confidence must be 
evaluated taking into account possible variations due to 
uncertainty in' the information available in the design. 

This v/ill help to measure the uniqueness of the design. 

Kov/ the problem of selecting a good design must be 
addressed. In the next chapter, measures v/ill be developed 
that will help to insure consistancy in our selection as 
v/ell as providing a reference from v/hich to base our deci- 
sions. 
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CHAPTER IV. 



SHIP DESIGN CRITERIA 



The existance of different commercial ship designs to 
accomplish similar tasks can be explained by the environ- 
ment and experience level of those involved in the design. 

The design selection process is greatly affected by these 
factors , because the real-world environment will influence 
individually each decision maker. Insight as to what deter- 
mines the merit of a ship design thus requires investigation 
of the decision maker's preferences. 

The objectives of the ship system which are to be 
satisfied by the decision maker's choice of alternatives 
are complex, however, a simple statement would list the 
commercial ship as primarily a means of improving the eco- 
nomic welfare of its owners by providing a marketable service. 

4, 1 Criteria 

The difference between the worth of ships is based on 
a number of criteria from v/hich the merit of a ship is deter- 
mined. Conceptually there is a single list of criteria whose 
various orderings reflect the preferences of the decision 
maker. The following figure attempts to give a sample cata- 
log of criteria. The determination of the relative impor- 
tance of these criteria is a function of the decision maker. 

34 



Table 1 



CATALOGUE OF CRITERIA 
Ref. [3,6.7,9,27,36,37,39] 

I . ECONOMIC 

Profits 
Durability 
Asset Life 
Obsolescence 
Salvage Value 

II. MISSION RELATED CAPABILITIES 
Transport Mission 

Speed 

Endurance, Range 

Capacity, Payload, Cargo Typ 

Patrol Mission 

Search 
Investigate 
Plant and Retrieve 
Secure 

Interaction Mission 

Engageinent 
Station-keeping 
Avoidance, Defense 

III. OVERALL SYSTEM CAPABILITIES 

Mobility 
Viability 
Flexibility 
Survivability 
Opera tability 
Ef f ectiveness 

IV. MARKETABLE QUALITIES 

Aesthetic Value, Style, Color 

Quality 

Prestige 

Personal Satisfaction 

V. RISKS AND UNCERTAINTIES 

Investment Risks 
Future Uncertainties 
Reliabilility 
Maintainability 
Availability 



The commercial ship ovmer, on one hand, is ultimately in- 
terested in economic criteria and in particular profits, 
whereas a warship's merit will be a strong function of the 
ship's mission capabilities. Where the ability of providing 
defense from attack may be important in the warship, this 
criteria is ignored in commercial ship design. Also, the 
economic criteria loses its eminence for a warship design 
because profits are no longer meaningful. 

Another difference between designs is reflected in the 
variation of risks involved in the design and the risk pre- 
ferences of the decision maker. The terra 'risk' refers to 
the possible outcomes associated with uncertain future events. 
For a commercial ship this may be in the form of unexpected 
breakdowns resulting in expenditure of funds and delays in 
operation, A risk associated with warship design may be 
represented by possible threats such as a surprise enemy 
attack. The risk preferences relate to the utility function 
of the decision maker. From this it is possible to determine 
the mechanism by v/hich the decision maker assigns weights to 
criteria, taking into account the values and risks involved 
for each alternative. The process of determining a decision 
maker's utility function uses successive lottery choices 
where the decision maker is asked to determine his preference 
between only two events. It is important to note that the 
decision maker is not always consistant in these choices. 
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The risks of the design are measured by many of the listed 
criteria. 

The multiattributed analysis of Keeney ^ and many 
others could be applicable in helping to quantify the owners 
preferences. For an example, if an owner of a commercial 
ship experiences a situation where a breakdown only occurred 
once during the life of the ship which required assistance 
to get into port, the risks may be acceptable. On the other 
hand, a single weapon impact on a warship may result in a 
total loss of the vessel and perhaps result in the loss of 
other units which it may have been defending. The risks are 
high and even though the event may be very unlikely, the 
resulting risk preferences of the decision maker may cause 
changes in the ordering of the design criteria, 

k,2 Measures of Effectiveness 

To be useful in a modeling context, each of the above 
criteria must have an associated measure. For those criteria 
which are basically quantitative this presents no problem. 

The economic criteria can be easily measured by dollars. The 
others, however, present a real problem, both to the designer 
and to the owner. This becomes most evident when it is nec- 
essary to make a decision based on tradeoffs between two 
qualitative criteria. It is not uncommon to assign measures 
to these qualifiers. The regulatory bodies have taken the 
liberty to associate floodable length with the safety of the 
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ship at sea. This is satisfactory if the proper response is 
obtained by use of such indirect measures. 

It is not the purpose of this paper to propose measures 
for the criteria. This is a function of the individual deci- 
sion maker. For the remainder of the paper, we will handle 
only the more accepted measures. We will use dollars as a 
measure of the economic criteria and as the measure of effec- 
tiveness for judging the ship design. Other quantitative 
criteria will be introduced as constraints, giving minimum 
standards for the design to satisfy, the freeboard requirement 
is such a constraint, 

4,3 Measure of Economic Criteria 

When making a capital investment decision it is impor- 
tant to identify all of the interactions that can affect the 
economic returns to the decision maker over the life of the 
decision. These interactions may be related indirectly as 
would be the case when the reputation of a ship affects 
the willingness of customers to do business. For this paper, 
we are restricting our interests to the direct economic 
effects, 

4,3.1. General 

In the context of the normative process, the decision 
maker attempts to tradeoff different criteria in terms of 
the possible changes in profit and other qualitative aspects. 
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Once the tradeoffs are made and the ovmer's requirements are 
identified, one economic measure of the goodness of a design 
is its ability to provide profit to the ovmer. In ship de- 
sign, this measure is subject to uncertainty because of the 
lack of information about revenues over the life of the 
ship. The measures in the literature assume that the reve- 
nues are unknown because a good forecast cannot be obtained, 
the result being that measures typically involve costs only. 
For the ship design, the costs can be predicted with a fair 
degree of certainty. 

If one assumes that the ship will be built and we are 
merely interested in determining the best designs out of 
the various alternatives, then the simplified problem can 
be analyzed by including revenues. A marketing study could 
provide an estimate of the income expected from the cargo 
transported. This would help indicate the relative earning 
ability of each design. The freight rate that would be 
obtained will be used to drive the solution technique in 
the sample problem. 

The o\vner's requirement for a desired cargo flov; rate 
would be calculated as shov/n in the follov/ing figure which 
depicts estimates of revenues and costs for shipping a 
different cargo transfer rates. The owner would select a 
value of the payload transported per year that would maxi- 
mize his profits. Along with this point estimate, he would 
also have an estimate of costs involved if the design 
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fluctuates from this estimate. This cost can be associated 
with a penalty cost to the designer for deviating from the 
owner's original requirement. Given the flow rate of cargo 
required, the o\vner is interested in the design that mini- 
mizes the cost of this service, subject to restrictions 
placed on the design. The owner would easily change his re- 
quirement on the flow rate of cargo if in so doing he would 
permit a design that could increase his profits. This gives 
the owner two options. First, he could direct the designer 
to design for a certain speed and cargo capacity, as is the 
present practice, or the designer could perform the evalua- 
tion to determine the new requirements. In this case the 
designer would need the criteria for setting up a decision 
model. This should require less work of the owner than to 
determine the owner’s requirements, and it would reduce the 
number of iterations of the proposed design process because 
the owner would not have to v/ait until the designer proposed 
each modification to the requirements. 

For the sample problem, the decision model is formulated 
in terms of economic criteria. By defining the general form 
of the revenue schedule, the cash flows can be determined 
by the designer. The model assumes a fixed freight rate for 
material transported betv/een tv/o specified locations. The 
designer would then be able to perform the tradeoffs in the 
design procedure concerned with the economic criteria once 
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the owner's preferences are identified. One way for the 
owner to project his desires is to identify the parameters 
for the calculation of the net present value of a design 
alternative, 

4,3.2, Net Present Value 

The Net Present Value (NPV) calculation takes a series 
of cash flows and determines an equivalent lump sum value at 
a specified point in time. It is based on the premise that 
there is a time value of money. The assumption is that a 
person will be indifferent to receiving one dollar today or 
one plus X dollars at some later time. This value is in- 
ferred from the fact that the one dollar can be invested to- 
day and have earned X additional dollars by the end of the 
period. We define the value of X as the discount factor. In 
general, the cash flows can be continuous or discrete, but 
for a continuous cash flow there will be an instantaneous 
discoimt rate. The cash flov/s may either be positive or 
negative. Also the discount factor may be a function of 
time. 

For ease of computation and also because of limitations 
on the projection of future cash flows, we will restrict the 
discussion to discrete cash flows and a constant discount 
rate. See Figure 6. Since data is available on an annual 
basis, we will use annual costs directly. The horizon of 
our calculations will be the expected lifetime of a ship. 
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It has been standard practice to use a twenty-five year life 
for ships. This has been substantiated by studies conducted 
by John McMullen Associated, Inc. for the U. S. Navy in 196?. 
At the end of the twenty-five years, the vessel will have a 
scrap value. This is usually a fraction of the initial cost 
and will be an input to the model. 

Looking at the operating costs when performing our 
present value calculation, we can make some important gener- 
alities. First, the operating level of the ship can be 
assumed constant for the duration of its life. This indicates 
that the crew hours and the fuel usage per year will be cons- 
tant. V/e can then calculate anniial operating costs if the 
wages and the price of fuel are known. If a current ma.rket 
price is assumed, then the discount factor could adjust for 
the effects of inflation. This may not be adequate if there 
is a restriction in the supply of fuel oil that causes it to 
increase in price faster than other expenses. One solution 
would be to utilize tv/o discount rates, one, for each different 
cash flow but for this problem we will assume only a single 
discount factor. 

Until now v/e have not included the maintenance and over- 
haul costs. These are significant when taken over the life 
of the ship. They should be estimated and included in the 
NPV calculation. These costs are dependent on the policies 
of the ovmer and data is not readily available. They will 
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be treated in this model as an annual cost whose value will 
be a fraction of the initial ship cost. 

The following list gives the major restrictions and 
limitations of our use of NFV, 

(1) Cash flows must be on an after tax basis, 

(2) Does not account for uncertainty in cash flows, 

(3) Assumes that there is no capital rationing, 

(4) Gives a poor measure for liquidity or timing 

„ 8 
preferences. 

Depending on the form of the capital rationing, modifi- 
cations to the straight NFV calculation can be made so that 
comparisons between designs are more valid. One method in- 
volves varying the discount factor and noting whether or not 
the relative merit of competing designs change. Another 
would impose a variable discount factor. This could model 
a temporary cut-off rate. For a strict capital rationing 
which limits the initial cost of the investment, the addi- 
tion of cost constraints to the design model would be another 
method of accounting for the rationing of capital. 

The other limitations of the net present value method 
indicates that it will not always serve as a good single 
measure of economic worth of the different design alterna- 
tives. In order to compare the results from this measure 
of effectiveness or objective with some others which are 
commonly used, the model calculates independently the values 
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of the capital recovery factor and the required freight rate 
for each alternative. This leaves open the option of selec- 
ting the objective to be used. 

4.3.3 Capital Recovery Factor 

The Capital Recovery Factor (See Figure 7 for Defini- 
tion) is the resultant of dividing the value of the annual 
cash flov/ by the initial investment. This can be used to 
find an implied or yield rate. If the cash flows and the 
investment base is the same as for the N.P.V. calculation, 
and the annual cash flows assumed are assumed constant, the 
results will be compatable with the N.P.V. measure. The 
capital recovery factor calculated is based on the total in- 
vestment and as such is not consistant v/ith the net present 
value measure. 

4.3.4 Required Freight Rate 

The Required Freight Rate (See Figure 7 for Definition) 
is that cost that must be charged to break even with the 
costs of operation, including depreciation or some return 
on the original investment. This calculation does not re- 
quire knowledge about revenues and is insensitive to change 
in tax rates. 
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= l.~ (2.-^3. -^5,) 

6 . 



RFR ( ) = 2. + 3, + 4. 

7 . 

Where : 

1. - Annual gross revenue 

2. - Annual operating cost 

3. - Annual loan payments 

4. - Annual depreciation allowance 
3. - Annual tax payment 

0. - Out of pocket investment 
7. - Payload weight in tons 
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CHAPTER V. 



ALTERNATIVES 



5. 1 Decision Variables 

The decision variables were selected to simplify the 
determination of the objective surface used by the optimi- 
zation technique. The model must incorporate variables that 
would first identify the physical form of the ship; second, 
identify the loaded condition, and finally identify the oper- 
ational capabilities. By using the linear dimensions and 
the coefficients of form, the geometry was determined. The 
specification of the draft and the displacement determined ^ 
full load condition. The only ship performance investigated 
was the speed of a transit. Thus the ship's speed was used. 
The following is a list of the decision variables, 

X^. ..... . .Displacement (A ) 

Xg. Length (L) 

X^. ..... . .Beam (B) 

X^ Draft (T) 

X^ Depth (D) 

X^ Prismatic 

Coefficient (CP) 

Xr, Speed (V) 

Table 2 

Is’ote that the midship coefficient is a dependent varia- 
ble. CM = CB/ CP. 
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5. 2 Design Alternatives 



The various alternatives that will be investigated can 
be identified as a set of the design variables. For each 
alternative, other parameters are held constant except when 
sensitivity to changes of these parameters were tested. 

This method of alternative identification was selected to 
maintain simplicity in the design model and to facilitate 
its understanding. Also areas of scarce information played 
a role in the final set of variables. 

The ship capabilities which have a major effect on the 
ship owner are the load capacity, speed, and endurance of 
the vessel. Prior to evaluating any design these must be 
known. Each contains a continuous range of possible values. 
For this problem, the endurance will be treated as a fixed 
quantity. The speed and load capacity of a given size and 
shape ship are now related. After the propiilsion plant con- 
figuration is determined, it is possible to evaluate the 
available space and weight for cargo. The size of the pro- 
pulsion plant is a strong function of speed for any given 
ship geometry. It should be obvious that the alternatives 
can be identified by the decision variables listed earlier. 

In the search for the better alternative designs, the 
feasible space must be identified. The following areas are 
discussed to identify the limitations of possible alternatives 
studied by the design model. 
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propellers and shafts in the design. In the model either 
a single or double screw arrangement may be selected. These 
restrictions are imposed by the program user. The model 
only adjusts the number of engines and the associated opera- 
ting levels and does not handle the actual engine arrange- 
ment, The differences in the v/eights of units other than 
the engines are not considered significant for the differ- 
ent configurations. There is, however, an adjustment made 
for the differences in appendage drag. 

5.2.2. Ballast 

As used in this paper, ballast is not considered part 
of the payload. An upper limit on payload for each alter- 
native would occur when no ballast is being transported by 
a proposed design. There may be a configuration where bal- 
last is transported. Of course, this is not an intuitively 
appealing situation, but there may be a configuration v/here 
ballast may be advantageous. This may occur if the ship is 
stability limited. If by carrying extra ballast an addi- 
tional container is transported, the benefits could out- 
v/eigh the costs. Even though ballast may be *off limits" 
to most designers and onerators, the only direct cost asso- 

X 

ciated with the ballast in a design with a fixed displace- 
ment is in the initial installation. The inclusion of bal- 
last impacts on the payload carried at the specified speed 
and endurance. For the purpose of this study, the cost of 

installed ballast will be taken as $1,000 per ton, 
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This 



is higher than the real costs experienced, thus shoiild help 

V 

to minimize ballast wherever possible, 

5.2.3 Full load Condition 

It is realized that designing to the fuel load condi- 
tion does not reflect all of the possible operating condi- 
tions v/hich exist. The program does not address directly 
the backhaul or other ballasted conditions. Through inter- 
action with the program, however, it may be possible to 
develop a design model that could incorporate such a ballast 
transit. This would require that the designer interact with 
the program and the results be carefully interpreted. The 
design model could in this v/ay serve as a small part of an 
overall system design optimization, 

5.2.4 Volume 

Using a similar argument as that for disregarding en- 
gine space location, the volume constraints of the problem 
were ignored. Except for the obvious case of determiping 
the container configuration, the volumes were not tested. 

A final hand calculation of volumes v/as accomplished for 
the final results. From these hands calculations, it was 
found that the designer may wish to be more careful in the 
handling of the volume constraints. First, it was thought 
the ship may be only v/eight limited for all but the container 
configuration. That v/ould mean that the space not used by 
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the container configuration would be ample for the volume 
of any other weight group such as fuel. 

5.2.5. Containers 

The variation in container sizes has been standardized 
since 1961. The American Standards Association standard 
container dimensions are given in Table 2 . In this study, 
only the 20' X 8* X 8' non-refrigerator containers were 
used. It was assumed that they would possess the strength 
necessary to stack the containers six high. This standard 
of the ASA is not accepted internationally, the International 
Organization for Standards (I.S.O.) calls for the stacking 

9 

of only four high. 

5.2.6. Terminal - Interaction 

The various interactions with the terminal were 
treated as parameters in the design. It was assumed that 
the ship relied on pierside handling and the average port 
delay time was not sensitive in the range of the loading 
capacity of the ship. This could have been modeled but 
again this level of detail was beyond the scope of the 
problem. There was one factor of the interaction that was 
incorporated. The model is constrained to a maximum number 
of containers per trip. This relates to the capacity of 
the terminal facilities. 
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5.2,7. Propeller Alternatives 



It was decided to handle single or double shafts se- 
parately. However, the effects on the overall propulsive 
efficiency due to different propeller designs were ignored. 
The problems associated with the propeller design could have 
been incorporated at the expense of simplicity. Instead, it 
was decided to test the sensitivity of the overall results 
to changes in an arbitrary propulsive coefficient. This 
would then indicate if more detail would be necessary. 
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CHAPTER VI. 
DECISION MODEL 



The objective of the optimization model will be to 
maximize the net present value of cash flows to the ship 
owners due to the initial cost and annual cash flows. 

Each alternative generated will be subject to restric- 
tions from rules imposed by the regulatory agencies, the 
restrictions due to navigational considerations, and restric- 
tions due to either data limitations or owners requirements. 
The following is a simplified mathematical model of the 
problem. The expressions are given in a simplified form. 

The decision variables v/ere covered in Chapter V. 

6.1 Objective function 

As determined earlier, the net present value measure 
for the capital decision problem could be one of the cri- 
teria used in making choices among alternative designs. Our 
decision model determined the net present value of the cash 
flov/s to the owner by calculating the initial building and 
subsequent annual costs. The objective of the program will 
be to maximize the net present value of the design. The 
objective function is not expressed solely as closed form 
equations, because some of the data is extracted from known 
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tables. 



6.2 Constraints 

During the design process, there are several con- 
straints that must be observed. These can be defined as 
either soft or hard constraints of the problem. They help 
determine the feasible set of design alternatives. The 
hard constraints refer to the constraints imposed by tech- 
nology and the laws of nature. These would include the maxi- 
mum draft limitation due to the channel depth and the limi- 
tations on geometry of the ship due to strength limitations. 
The exact value of the hard constraints are sometimes diffi- 
cult to determine. For example, the limit of strength is 
not well defined. In these cases, rules of thvunb have been 
developed. These are the result of attempts by the regula- 
tory agencies to quantify the factors involved. 

The soft constraints, on the other hand, have no phy- 
sical significance. They result, not from ship design 
alternative, but from requirements of the owner or designer. 
These could be imposed in order to insure compatibility with 
other elements in the total design of the ocean transporta- 
tion system. The maximum constraint on containers permitted 
per trip is an example of such a soft constraint. It could 
also be a limit indicating the range of data or experience. 
The limit is obvious in the powering calculation. The soft 
constraints thus limits the set of feasible design 
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alternatives. Even thoiigh there may be good alternatives 
outside the experience of a designer he may not have con- 
fidence in its potential. For this problem, we have im- 
posed soft constraints due to the availability of data. 

If these become tight constraints during the final itera- 
tions of the design then a direction will at least be iden- 
tified for further investigations. 

For this problem, we have segregated the constraints 
into four areas t Those due to the requirements of regula- 
tory bodies, navigational restrictions, interaction restric- 
tions, and data restrictions. (See Figure 8) 
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MATHEMATICAL STATEMENT OF THE 
CONTAINSRSHIP DESIGN PROBLEM 
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CHAPTER VII. 



OPTIMIZATION TECHNIQUES 

The optimization problem as stated in the last chapter 
is a mixed integer, non-linear, mathematical problem with 
some expressions not given in closed form. The relation- 
ships used to determine the powering requirements of the 
ship are a combination of working graphs and tables. The 
relationships betv/een several decision variables are in- 
volved and any simplifying assumption necessary to obtain 
a closed form expression would introduce unacceptable errors. 
Without such a closed form solution, analytical solution 
techniques cannot be incorporated. For a survey of the 
approaches applied to this problem, see Appendix to 
Reference 51. 

One method of determining the powering requirements for 
design alternatives has been tabulated in working graphs. 

The Taylor Standard Series of model tests is an example which 
presents the powering requirement as a function of five varia- 
bles. (Displacement, beam to draft ratio, volumetric co- 
efficient, and speed to the square root of length ratio). 

There have been attempts at determining the analytic form 
of the relationship, hov/ever, one or more of the independent 
variables are usually set equal to constant values. These 
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variables selected are indicative of a special type ship. Be- 
cause the relationships are not available for containerships 
in general, these analytic forms were viewed as being too 
restrictive. The solution methodology used must be compati- 
ble with the data in its available format of graphs and ta- 
bles, The identification of an objective surface is possible 
using the available information. Also the constraints can 
be readily evaluated. This situation suggests that some form 
of search routine will be able to locate the location of the 
optimal objective value which meets the constraints. 

7.1 BASIC SEARCH ROUTINES 

7.1.1 Random Search 

The Random Seach Technique applied to the ship design 
process is covered by I>iandel and Leopold The method in- 
volves searching each variable sequentially, each time sam- 
pling at random from the acceptable range of the variable. 

As the procedure progresses, the search is concentrated 
around the more attractive solutions. This results in the 
determination of a solution for the first stage. The pro- 
cess is iterative, it is started over many times, with the 
best solution of previous iterations being retained. How- 
ever, there are drav/backs to its use. First, there is 
the large amount of computational effort required. Secondly, 
there is difficulty in identifying the feasible range for 
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each of the variables searched. It maps out the surface dur 
ing the process which provides additional information to the 
designer. Unfortunately, its randomness makes it difficult 
for the designer to interact with the program. The random 
nature does, however, help to insure adequate coverage given 
enough iterations. 

7.1.2 Climbing: Techniques 

The Hill Climbing Techniques employ the information 
gained from the gradient of the surface to determine the di- 
rection for the next move. This technique is often used in 
conjxmction with other techniques which determine the dis- 
tance to jiinp. The methods require an initial point and the 
derivatives at that point. The process is usually slow to 
converge due to possibility of rapid oscillation in the move 
ment vector as the solution is approached. Problems due to 
special surface forms have been handled by variations of the 
technique. Many of these include systematic rotation of the 
coordinate axis.^^ The laajor difficulty v/ith this technique 
for our application is the problem of keeping the search 
within the feasible region. 

One method of insuring that the search remains in the 
feasible region is to impose severe penalties for violating 
the constraints. The Sequential Unconstrained method uses 
this barrier technique successfully as discussed in the next 
section. 
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7.2 SEQUENTIAL UNCONSTRAINED MINI^:IZATION TECHNIQUE (SUMT) 

Ref [15] 

7.2.1, General 

The name of this technique is very descriptive. It is a 
process which transforms the constrained problem into a se- 
quence of more easily managed problems without constraints. 
The set of solutions to these problems possibly converge to 
the solution of the original, constrained problem. There is 
no single method of determining the sub-problems nor the 
associated solution method. Presented here is one method 
that has proven useful. 

The SUI>IT is designed to solve the optimization problemi 




SuBrecT to: >0 .1^ 

WueRF X • • • j Xk)"^ is ^ n-bine-osioA; colum/o vfcto/?. 

The technique applies a minimization method to a modified 
objective function P(x,r), where r is a given parameter, 

= f(j?) /r 

From the function P(x,r) one can associate the additional 
terms to a penalty. Depending on the value of "r", the pen- 
alty takes different forms as illustrated in Figure 8 . 

In the limit as "r" approaches zero, both penalty terms 
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represent barriers. The systematic reduction of the value 
of the parameter "r" and the solution of the associated 
modified objective functions, under suitable conditions, 
will result in a sequence of minimums approaching the 
solution of the original problem. 

P(x,r) = X'Cr) 

HM. X'(r) = 

Lim. fix'cr)]- f ^ V' 

V* -V o 

VVHt=Rt^ >>' - TmF SOi-UTtOA; THr 

1 2 

The suitable conditions are as follov/s: 

(1) The feasible region is non-empty. 

( 2 ) f(*) amt. 1,2^... m are convex and 

continuous and , bi-vz, - • m+r 

are linear function. 

C3) For some constant K » "the set 

is bounded, (This prevents a "minimum 
at infinity " ) 

( 4 ) P(y^r) is strictly convex for all r>o . 
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The objective surface of interest is not necessarily 
convex. This means that the results of the technique may 
be local optimums as opposed to a global solution. It must 
be emphasized that there is no guarantee of finding the 
problem solution in a single iteration. Additional runs 
with various starting positions would be necessary so that 
the characteristics of the surface could be fully explored. 
Only then would confidence in a single solution be justified. 
In this thesis, the application of the technique to a non- 
convex surface will be successfully demonstrated. 

There are a variety of ways to implement the SUI'.iT 
method. The one utilized here can be described as follows. 
Given a starting point, the method determines if the point 
is feasible. If not, the technique is applied to find such 
a point by using an objective function in the first phase 
v/hich measures the degree of non-feasibility. After a 
feasible point is found, the value of the modified objec- 
tive function F(x,r) the first, and the second derivatives 
with respect to the decision variables are determined at 
that point. The direction for the next move is calculated 
using the Newton-Raphson r.:ethod. If it is not possible to 
determine the inverse of the Hessian Hatrix (See I-age 69), 
the matrix of second partials is perturbed, and a move is 
determined when the matrix inverse is able to be calculated. 
This results in an "Crthoginal Hove" as identified in the 
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program. Once given the direction of the move, the Golden 
Section Search is employed to find an improved point in the 
move direction. This process is continued until the mini- 
mum of the program is found. This completes the first sub- 
program, Next the value of "r" is reduced and the process 
is repeated. This continues until the sequence of solutions 
converges. The result is the solution of the original prob- 
lem with constraints. 

In the next sections, the Newton-Raphson Method and 
Golden Section Search is discussed in more detail. 



7.2.2 Newton Ranhson Method Ref [52] 

This method is used to identify the move vector. It is 
an indirect method for solving simultaneous non-linear equa- 
tions, It utilizes the necessary conditbn for a point x 
to minimize a function Y(x). All the partial derivatives 
must disappear at that point. 
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The Nev/ton-Raphson Method then solves -the system of 
equations by using the derivative to estimate the functional 
value in the neighborhood of a point. This is done by per- 
forming for each equation a Taylor Series expansion of 

about X . If the higher order terms are ignored, 
the following equations are obtained, 

Yi - Yj (■?*') +Z Y-p j j IT 

J P=. 

X - ( Xi,v y ) ^3,k j Y ) • • ■ 
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The determination of requires the solution of 

these simultaneous equations. If vie define the Hessian Ma 
trix of second partial derivatives j 
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For the determination of the move vector, we are inter- 
ested in the direction along which to search for a better 
solution. The direction vector is defined by (vY 

for each step. The next step involves calculating the dis- 
tance to move. The procedure is then repeated until the 
functional values of become small. This process may 

not converge to a unique solution as seen by the graphic 
example given by V/ilde and Beightler ^3 in Figure 10. Assu- 
ming the system of equations has only one member, then this 
curve could represent the functional value of an equation. 
The roots or values of the independent variable X for which 
the function goes to zero is represented by the intersection 
of the curve and the X axis. In this case, there would be 
tv/o possible solutions at points a & c. Depending on the 
starting position for the Newton-Raphson Method, the con- 
vergence to a root (solution) takes on different character- 
istics. The zones of definite, possible, or non-convergence 
are identified. The program using this technique identi- 
fies the non-convergence condition v/hen it occurs, 

7.2.3 Golden Section Search Ref [52] 

The Golden Section Search is a method for determining 
the optimum of an objective function by successfully elimi- 
nating regions for investigation. This method is used to 
search along the direction determined by the Newton-Raphson 
.Method to find the best solution. This new solution 
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Figure 10 NEV/TON-RAPHSON PROCEDURE 
Ref [ 52 ] 



point is then used as the starting point for the next itera 
tion in the minimization of the objective. This method was 
originally designed to search a set interval. The SUIvlT pro 
gram determines this interval in each case by searching 
along the vector until a functional value larger than the 
initial one is found. This defines the search interval. 

The upper bound is determined by looking at successive 

points imtil a larger value than the initial one is foxmd. 
The distance between successive points is taken as Ou. . 

0; = Z 

i-o 

then 

(x, -e ejs)> /(X.) 

where 0^ is defined with the smallest non- 
negative integer which satisfies the above 
equation and S is the unit vector in the 
move direction. 
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After the interval is defined, the search proceeds to 
determine the optimxim value. This is accomplished by divi- 
ding the intervals into two parts and looking at the deri- 
vative at the intermediate position. This would indicate 
in which sectio.n the optimum v/ould most likely occur. The 
process v/ould continue until the optimum is foimd. The 
method which minimizes the number of iterations without 
additional information can be shown to be a method of 
partition which divides the remaining interval into sections 
of the ratio T (1 . 618033989) . This technique converges to 
a unique solution provided the objective function is strictly 
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unimodal as defined by Wilde and Beightler, 
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7.3 SUWT COr/IPUTER PROGRAM 



The optimization technique previously described has 
been programmed by the Research Analysis Corporation as out- 
lined in Reference l6 . This program requires that the user 
provide four subroutines. These include first a subroutine 
that reads in values as needed by the user in the other 
supplied subroutines; secondly, a subroutine which defines 
the value of the objective function and constraint func- 
tions; thirdly, a subroutine v/hich determines the partial 
derivatives, and finally a subroutine which provides the 
second partial derivatives of the objective and constant 
functions. 

The program was vnritten with the capability of linear 
interpolating for the required first and second derivatives. 
This required that the objective function be developed from 
the design model. Such a routine would provide the required 
objective value for any value of the input variables. This 
routine thus describes the objective surface upon which the 
optimization technique v/ould search and is used in the in- 
terpolation of derivatives. 

The development of any objective routine must insure 
that for each combination of the input variables, there is 
defined only a single objective value. The program used 
defines several surfaces, one for each set of input para- 
' meters. The tabular .inputs for powering calculations v/ere 
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easily handled in this method. Simplified expressions for 
these would lead to the use of closed form expressions which 
could then "be handled in other optimization techniques, one 
of w’hich would be geometric programming as used in Refer- 
ence 10, The model developed in the next chapter attempts 
to incorporate sufficient detail during a preliminary de- 
sign stage so that the full capability of the SUI.IT optimiza- 
tion technique can be explored. 

The program requires the selection of an initial 
starting point and various control parameters. The author's 
experience v/ith these and other features are covered in the 
concluding chapter. The reader is directed to the user's 
manual, Reference 4-3. A complete explanation of ths appli- 
cation of SUI.iT v/ith input, set-up, and a sample run is 
covered in the Appendix, 
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CHAPTER VIII. 

DEFINITION OF THE OBJECTIVE SURFACE 



8. 1 General 

The strength of the search technique depends upon the 
explicit nature of the objective surface. The selection of 
the decision variables and design relationships has focused 
on the need to determine an explicit, single valued objec- 
tive as a function of the decision variables. The general 
application of an optimization technique as described ear- 
lier in this paper requires the determination of such a sur- 
face, The actual surface defined is intended to be repre- 
sentative of an actual decision, hov/ever, the main interest 
here is in obtaining a model with which to test the optimi- 
zation technique. 

In this chapter, the containership design model which 
is used to determine the objective surface will be described. 
The actual relationships and other data used has been com- 
piled in the Appendix. The methodology used to construct 
this model may assist in the understanding of its final form. 

The first step involved defining a tentative set of 

variables. The decision variables introduced in Chapter V 

represent the primary, independent variables chosen. It 

would be possible to define the surface in terms of only 

these variables, but the comnlexity of the expressions would 
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lend the model awkward. Instead of following the usual pro- 
cedure of introducing these new variables v/hich are depen- 
dent upon the primary variables and then structuring the 
objective function in terms of these secondary variables, a 
more convenient method was used. This also helped to mini- 
mize the number of secondary variables defined. The trans- 
formation defining secondary variables were viewed as "black 
boxes" or subproblems. Often explicit relationships between 
the output and inputs \vere available. At other times, the 
relationships were given in terms of tables or graphs. 

The construction of the model started at the end where 
the output objective v/as desired. Using the subproblem 
technique, the model was built up from the objective end. 
This systems approach facilitated the model construction. 

As subproblems v;ere defined, their inputs became outputs of 
the next generation of subproblems. The goal in defining 
the subprograms vras to determine the inputs necessary and 
the relationships required to obtain a simple and single 
valued relationship with the required output. It was found 
efficient to make extra subproblems v/henever the relation- 
ships betv/een the variables could not be expressed in a 
simple, single, closed form algebraic expression. This 
v/ould provide the definition of new variables. Another ad- 
vantage of proceeding in this manner v;ould be the capability 
of having several groups working simultaneously on the 
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model after the first few steps. The assignment of subprob 
lems should prove to be manageable in a real design effort. 
The following black box exemplifies the subproblem tech- 
nique , 



The process culminates when inputs to all the subroutines 
are either the primary decision variables or outputs from 
other subroutine. Only as the process is near completion can 
the final set of decision variables be identified. 

The following block diagram depicts a single design iter- 
ation. The diagram represents only a small part of the total 
computerized program. It shov/s the development of the ob- 
jective surface in terms of the design variables. The output 
of the program is then used in conjunction with the SUIuT to 
obtain new solution points as shown diagramaticly by the feed- 
back loop. Evai -though the Net Present Value objective was 
selected, it would be a simple matter to change to some other 
objective for the problem. See Appendix A for definitions. 

The relationships that are used in the black boxes can be 
found in Appendix E by their assigned number. The program was 
designed so that the user could substitute other relationships. 
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INTERACTION BETWEEN SUMT 
AND THE SHIP DESIGN MODEL 

Figure 16 



I 





There are three blocks in the flow graph that require exten- 
sive calculations and additional explanation. These are the 
blocks which determine the required propulsive power, the 
power plant selection, and the container configuration, 

8.2 Required Propulsive Power 

Given a specific ship's form, methods have been devel 
oped to predict the required propulsive power to drive the 
vessel at different speeds. The process involves the use of 
model tests and correlations involving past design experience 
The best approximation requires that a separate model be 
tested for each design alternative investigated. This is too 
expensive a procedure to follow when there are several differ 
ent alternative geometries. Over the years, systematic stu- 
dies of the resistance properties of ships have been conduc- 
ted, These have resulted in several hull forms being tested 
and the results compiled in working graphs. Each of these 
series of tests typically limited the number of independent 
variables by fixing the value of one or more shape coeffi- 
cients. It should be obvious that the use of any one series 
to determine the power requirement for an alternative would 
be only an approximation if the geometry of the design does 
not match that of the model series. 

In the program developed, the Taylor Standard Series 
was used. In this series, the relationship between two of 
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the decision variables is fixed. For the Taylor Standard 
Series, the block coefficient is related to the prismatic 
coefficient as follows: 

It was found necessary to impose this as a restriction in 
the mathematical model because of the strong relationship 
betv/een the block form and the revenues generated by the 
design. When this constraint was not imposed, the program 
would select a box shaped midship section to house the con- 
tainers, The resulting changes in pov;ering requirements 
could not be calculated using the series data. Even if ano- 
ther series were used instead of the Taylor Standard Series, 
the program would experience similar difficulties with other 
pairs of decision variables. This constraint reduced the 
number of decision variables by one, 

8, 3 Discrete Power Plant Selection 

This model assumes that given a required SHP, a pre- 
ferred engine configuration can be identified. The problem 
of discrete power plants experienced with gas turbines has 
been incorporated into the program. Because the speed and 
distance values remain constant for the selection, the gas 
turbine engine alternatives are ordered according to their 
initial cost and speci‘’''ic fuel consumption. An engine con- 
figuration is identified by the number of gas turbines in- 
stalled, In this problem the GE built L.M. 2500 is used. 

As the SHF increases, there may be discontinuities 
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in costs. Because the gas turbine is available only in dis- 
crete sizes, there exists many situations in which the in- 
stalled turbines may not be operating at maximum rated power. 
The following decision rule is proposed for the selection of 
the number of gas turbines for installation. 

It is assumed that the initial cost of the power plant 
monotonically increases with the SHP required (See Figure 21 ), 
The variables cost for fuel varies proportionally to SFC for 
a fixed required SHP (See Figure 22), Because of the nature 
of gas turbines, their SFC decreases monotonically as the 
pov/er output is increased. It is, therefore, desirable to 
install turbines that operate near full power. It is assumed 
for this model that the best engine configuration is the one 
that is just able to meet the SHP requirement. Any combina- 
tion calling for fewer engines cannot provide the power re- 
quired and is thus non-f easible, A combination that operates 
extra eng.ines v/ill run at lighter loads and thus at a higher 
SFC. The desire of minimizing initial fuel costs associated 
with the power plant selection is consistant with this deci- 
sion rule. This ignores the change in maintenance costs with 
the variation in engine loading. 

Past history has demonstrated the reliability of gas 
turbines and it compares with other types of propulsion units. 
Typically more units have been designed into the plant as 
reserve pov/er sources. This increases the reliability of the 
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total system at an increase in cost. This model permits the 
designer to require additional standby units as desired, 

8,4 Payload Configuration 

There are two steps in determining the configuration of 
the containers carried on board a ship. The first determines 
the available space and v/eight available for the payload. The 
second determines the loading so as to meet the weight and 
stability requirements. This reduced the number of constraint 
programs. 

The space is determined by obtaining the maximum dimen- 
sions of the available space for the containers, then divi- 
ding the ships into decks. The weight available for the pay- 
load is determined by incorporating Archimedes Principle, The 
payload and ballast weights must equal the difference between 
the displacement weight and the sum of the other weights pre- 
viously determined. The program then loads each deck sequen- 
tially until one of the constraints is met. The number of 
containers per deck is calculated as follows. The containers 
carried on decks below the v/aterline were treated together and 
assumed to have a shape factor relating directly to the ship 
block coefficient. For those containers carried above the 
waterline and above the main deck, each level was considered 
to have a shape factor equal to the coefficient of the water 
pXane area. The program assumes that the coefficient of the 
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water plane area can be expressed as a function of the pris- 
matic coefficient. The shape factor is used to determine 
space for containers on each level by multiplying the shape 
factor by the dimensions of tlie cargo space. 

The model test for the stability constraint uses a cal- 
culation of the available v/eighted moment at the designed 
displacement . The program starts by placing all of the 
available weight as ballast in the double bottom of the ship 
without regard to space limitations. This automatically in- 
sures that the weight requirement is met. If this results 
in a feasible solution, that is, if the non-linear stability 
constraint is not violated, then small increments of weight 
are removed from the ballast location and placed in the form 
of containers in the cargo hold. The only change in the 
weighted moment v/ould then be due to the movement of the 
v/eight in the form of ballast in the double bottom to a 
position at some higher position that v/ould correspond to 
the position of a container. This process is continued se- 
quentially until the actual full load GM is reduced to 5 per 
cent of the beam or one of the restrictions is encountered 
resulting from stacking the containers. ?or ease of calcu- 
lation, the model groups the containers. Each deck is filled 
one tenth at a time. The program starts with the lov/erst 
possible level and fills each deck sequentially. The model 
v/ill stack containers external to the main deck. The 



86 



designer can limit the number of levels in this configura- 
tion. The sample problem assumes that the maximum number of 
containers that can be stacked is six. This is consistant 
with the ASA requirements for this size container. Also, 
the number of levels above the main deck v.'as limited to two. 

8.5 Input Parameters 

Associated with many of the subproblems shovm in the 
flow chart are parameters that must be provided to complete 
the calculations. These constants may be supplied by the 
program user. The input parameters are listed in Table 4. 
The adjustment of these parameters permit the designer to 
approach the real design problem and to introduce the 
effects of other criteria expressed by the ship owner. 

The ship availability refers to the per cent of time 
that the ship will be operational. The design and cost 
weighting factors are provided to permit the user to modify 
the equations used in the computer program. 
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Table 4 



PROGRAM INPUT PARAMETERS 



PARAMETER 


SYMBOL 


VALUE 


DESt&M ?/«\K/^MGTtRS 


, 




PROPU COeFFidlGAJT 


PC 




GK)l>UT{AKfCt 


£ 




DiST/)aJC6 PORT^ 


DBF 




SH(P /\V/)tL/lBl(-(TV 


AV-4IL 




hJUltA)i&A or H 4 n S 


NSHA ft 




NO. OF B\)(TiAje‘^'^ 


N£f<Jeu 




wo. OF SP/'RF e/W6-lA;FS 


A/eXTRA 




HA)<. rl4'^(.c>4b WT. 


WW 




M»/u. PAVlo41> V/T. 


WPMIAJ 




Max. <ioa/tai/v^4 stack(>v6- 


fJV 




MAX. STACKI/VO ABoi/t DK. 


MX) 


• 


AVFf^AC>F CO/J7A;vFN. '\/r. 


WC 




Bl=FecTlO& OoajT. LBAJCFUi 


BLENCrT 




BTtF<n:\WG CoxfT- uyi'DTH 


£wi*dth 




S rFECn OF Q O^JT , T)eTTf/ 


Et)FPTN 




‘DouBlG' Wrro/^ Heienr 


'DoUBL.F 




B/lLLAbr We[(rH^T KO. 


K^rB 




Max. 5ffp ro/<’ 


SI4PM 




$F/q SPffT) FActoa' 


BftR L 




5 PC. PAC.TC3/^ 


BetAE. 




Ml SC. W/Ft(rHT 


wxx 




MlSC.. uypiG-HT KG-. 


XK6-X 










OfSKoAJ u/to(^^(4T^A;6- F 40 TO/^S 






p4C.To/^ OwTF rr 


- v/Wo 




FACTOR FOR stffc wr. 


wws 




Factor for Nvku-HA/eR'fwj 


wwm 





* tJe6-A,T(0<F VAtwe (A;n><c4re5 p,qoPULUoAj , 
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Program 



Table 5 
INPUT 



PARAMETERS 



PARAMETER . 


SYMBOL 


VALUE 


NAVICATIOM CoW5TRAl*/r5 






MAX. LEM6TH 


LEn&TM 




MAX. b£AM 


beamm 




MAX. D«AFT 


DRfl FTM 










economic . PA'MtAFfeRS 






FRFI6-I4T RATF 


FR 




QALL F££ 


PCF 




P/HLY 'PoU.T F££ 


V>P F 


• 


PRiCt OF Fuel 


PRICE- 




AhJkJUAu. IajSURAVCC f^A-TE 


AIR 




ASSET C(Ffc 


NN 




SALVAGh VALUE 


SV 




■ discount factor 


D.F 




■ Subs 'bv Rate 


SUSIDR 




PER C.F a/7 P' ( N A/V CE'P 


pf 




IRTEREST F<AT? 


RAre 




Coat, tax rate 


CTR. 




INVF3TMEAT T4X C/?tDlT 


AfTC . 










Cost vveiehT/a/e factors 






FACTO/^ PO/i OJTFCf CO^T 


CCo 




F/hsto^ fo^ Cov> ( 


CCS 




FA<J0A FOA. MAUH/je/C'{ c. 












PROGRAM PAAAMETCRS 






OBjccrit/c SEc£k-t(o^ 






CC^STAAi^T 5£K/SiT(i/lT'( 


X>eLrAr(X) 
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CHAPTER IX 



ANALYSIS 



Responding to the needs of the owner as outlined earlier 
the optimization techniques are applied to the ship design 
model. The results of the investigation are presented to de- 
monstrate the capabilities of the computer programs. Prior 
to discussing the actual problem solution, an attempt was 
made to validate the design program by comparing the program 
results with three containership designs. 

9. 1 Program Validation 

The design program was used to develop details of three 
separate preliminary containership designs. The first design 
was developed by George G. Sharp, Company, during the evalua- 
tion of gas turbine"^power plants. The remaining two designs 
are products of the Maritime Administration CMX Project. The 
designs were developed separately by Bath Iron Works and the 
Newport News Shipbuilding and Dry Dock Company. The follow- 
ing tables indicate the similarity in results. In each case 
the design program v/as given the principle dimensions of the 
ship. The input parameters used in each of the designs corres 
ponded to the values obtained from the respective preliminary 
designs. 
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Table 6 

MODEL VALIDATION 



V/lRiABLes 


G.& Sharp 


Model 


bisPL/ictMewT 


'50i,5-<\ 




Le/NltTH 






beam 


^3 


^?3 


DRA Fr 






DfPth 


FFo 


rs’. o 


P'^ISFIAT/C C-G^F. 




<65' 


15Loc_k Coff. 


-Ffc 1 


,60 


Volume Cofp. 


• oo;^r 


.oo;?^ 


Vf L.OdTY 


S’. 8 




Rj^l wt. 


lA >7 0 


J03S" 


MA-CH. wr. 


112 5 


/f 


OCAf-FfT WT- 


( fer'i' 


;f6o 


STFFl WT. 


9oX0 


gT/A 


Misc. wn 


X<o% 


3oo 


Ballast wt. 


o 


/;? 8fc 


T^^Yi-oAb WT, 


1 S’ 3 V S" 


I'if.sn 


64A 


;.T3 


9. 1 


No. <l 0 WTAIAJ£(«S 


75*/ 




To'Pji't.t 




3 


BEUOv/ t>K«< 


wA 




No. Sfaftjs 






TPY _ _ 


WA 


n.s" 


ShP Rfq 


5SFOO 


S' S' 6HO 


CdstI 


FS’.FS 




Aod. 


2,5' ^ 


AM*f 


RFR 


NA 


..U7.S" 


C4RF 


MA 


. !To 


FR 


W A 


S"c?. 


NW 


^)A 


I'^.SS 




i 


i 



«e^L?cr« «•<=- '*^ 



Table 6 



• MODEL VALIDATION 



V/ARlABLFS 


NFvVPo/2 1 
NBWi 


Mot>ec 


B/^TH 
1 /5o// 
Wo/^/c'5 


HoT>5<- 


DlSFLACeMENn- 


hXloo 


34400 


51 l7o 


3(t^o 


LENetH 


696 


69t 


6AF 


64s- 


Beam 


lo3 


los 


86 


86 


DR4fT 


RT.S" 


49, S 


31.6 


31,5- 


DETmi 


6o 


6o 


FA. 8 


5^.7 


Prismatic- c- 




■S'-ts 


.639 


•439 


ftLOEk COEF. 


• S33 




.649 


.629 


Yolum. CoSF. 


.0033 


.0033 


.009 


.ca^ 


VEUo C-ITY 


43. J 




40.9 


ZO.H 


Fufl wr. 


32oo. 




36/ O 


inie 


K/VCH. v/n 


iU2. 


1282 


A/4 


1 o73^ 


OUTFrr WT. 




/6 30 


MA 


36oZ 


STFtL WT- 


Rixs, 


89 AT 


A/A 


ies9 


h i^.C,. WT. 




“3)00 


NA 


Zoo 


i54LL4srwr. 


367. 


IT90 


'"5ioo ^ 


S' ^>oo 


PaYeoAT \VT. 


I(^ Okv 


15-038 


NA 


fo^lo 


6M 




6.43' 


NA 


^,37 


N/o. CoAJTAive/Ci 


1 ^ s~o 


1268 


^2? 


943 


roTsit>f 


X 


3 


2 


3 


BicCcw *L>cCf<: 


7 


7 


G 


6 


s/c. SmA(T3. 


( 


( 


1 


/ 


TPY 


fJA 


(S-.g^f 


KlA 




Sa-p ffs. 


Zyooc) ! bf^oo 


.^Sooo 


ZA loo 


Co^Tl 




18.1 




22.9 


4 oc 


UA 


4,1 


wA 


4.F 


RPR n 


/M 


;ii4 


AJA 


37.93 


CRP 


A/4 


.5-t 


lUA 


,/4 


FR 


RA 


?o. 


NA 


Fo , 


AjpV 


RA 


17.6 


NA 


4.9 




Z 


£ * 


3 


3 



^ /2>V sr^ociiATC. 

^ ^/H-A/V^T LJ trfC^H’V rO^ SHO/^T<r^ l/6T^S<0-x4 
fj^ = aJoT 
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9. 2 Program Inputs 



The application of the decision model to the scenario 
requires several input parameters. Unless otherwise stated, 
these parameters were used for each of the designs analized. 
The parameters are categorized as either design, economic 
or financial. The design parameters include the assumptions 
necessary to determine the ship design as well as several of 
the constraint values. The economic parameters cover the 
various costs and revenues generated by the design. The fi- 
nancial parameters relate to other economic aspects that in- 
directly relate to cost by indicating the owner's preferences. 
Table 18 gives the values of the coefficients used in the sam- 
ple problem. 

9. 3 Procedure 

The application of the minimization technique requires 
the identification of a starting position from which to 
search the objective surface. The primary concern is to start 
each search v/ith a feasible solution to the problem. From 
this point, the SUIVIT Program will take over the procedure and 
then locate a local optimum. Because the objective surface 
is not always convex over the feasible domain, other starting 
points or solutions may result in the search culminating at 
different local optimums. The concern is then to start the 
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program a sufficient number of times, at different locations 
so that the surface characteristics can be deduced and the 
best local optimum can be taken as a solution to the prob- 
lem at hand. 

The procedure was to select four starting solutions, 
each different, so that the total search would encompass 
most of the feasible surface. The first starting point, 
which is labeled Alternative 2 on the next page, represents 
a modification of a design used by George G. Sharp Company. 
This alternative has similar characteristics of the actual 
design which is -shown as Alternative 1. The major differ- 
ence is in the number of shafts in the design. 

The second starting solution was chosen with a higher 
displacement. Alternative 3 was used for a starting point 
even though it had been an output of a previous optimization. 
The third and fourth starting points, double and single 
shaft arrangement for a smaller ship. They are shown as 
Alternatives 5 and 7. 

The remaining alternatives show the results obtained. 

Not shov/n are two attempts to start the program at non- 
feasible initial positions on the surface. This resulted 
in an unsuccessful attempt by the program to find a feasi- 
ble starting point, ^or the nonfeasible alternatives, the 
measure of effectiveness v/as most easily minimized by re- 
ducing the size o-^ the design. No revenues v/ere generated to 
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IDENTIFICATION OF ALTERNATIVES 



\ 

\ 

\ 

VARlABLES'^^ 




© 


© 


© 


DI5PLACEHENT 


30,«>5‘?. 


30.659. 


37,6mc 


H53 io . 


Le^76TH 




750 . 




7fei. 


BEAM 


^3. 


^13, 


7m. f 


A8.S 


DR/\FT 


3'7.3 


;^^.3 




33.9 


DETn-H 


55. 


55. 


5"g. 7 


5*7.4 


P/Rl5N1ATnC Co. 


.65' 




.6-7 


.47 


Block co&f. 


.56 


.St 


.43 


.t3 


Youum. CoeF, 


.OOZ.G 


. 007.5 


.oo3t 


.oo 34 


V'Fi,OCITV 




25. T 


22.8 




Fueu WT 


ttAt.3 




5*7feoa 


^ O 13.7 


K4cH. WT 


Io7o.3 


1 ( H7.3 


76M.7 


1 12?.’? 


OUTFIT V/T. 


?33o.3 


S53o.^ 


S'! 7.i .7 


1 


STEEL WT. 


snis-.<i 


S7i5,'^ 


8U6.3 


?r?3, 


>ifSC. V/T. 


l>oo. 


3oo. 


3oo. 


3oo 


S/?LL/^5T WT. 


aj.6 




^65.7 


^o^M,9 


P/^YLOAT^ VVT 


y.STJ.A 


7,? 7 a % 


1 S,^-tZ. 


1 




■ (5-.S 


/6.f 




8.07 


Ko. Co/srrAi^ms 


5-I8 


475 


/ (25 


/ / 


roPs<X>e 


O 


o 


Z 




Btcovo 


vT 




4 


4 


Wo. .5H4FT5 


i 


a 


2 


2 


KO. BjCtIWS 


3 


4 


2. 


4 


SH? RfQ. 


57, 2JI. 


4l, 63S. 


7F.J77- 


54, 5°o 


C05T i 


3I.Z6 


3 2.U, 


F7.38 


■ 


AoC 






7.78 


/Z, M 


/^FR 


82.37 


/oo. tr 


38.33 


Y3.T0 


CP.F 


.13 




.^S 


.57 


FR 


fOO. 


|0O. 


(oo . 


/oo. 


MPY 


Io.h7 


o*~ 


48.08 


S-/.70 


Pro&rak 


NA. 


IV Put 


OUTPUT //A/Pur 


ouTptrr 



4 t*/l>JCATe5 NACkTu/f 

Table 7 
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IDENTIFICATION OF ALTERNATIVES 



"-.^creRMATivE 

\ 

\ 

\ 

VARIABLEi'v^ 


©, 


® 


© 


® 


DI5PL/5CEMCWT 


IS,ooo. 


Zl.llt. 


if^OCO. 




length 


Goo. 




coo. 


67o, 


BEAM 




'78.2 


50. 


T?.Z 


DRA FT 


zsr.o 


18. H 


2f.o 


32.0 


DEPTH 


AS. 


AS. 


^5“. 




PRISMATIC- to. 


. bo 


.bo 




.62 


Block coef. 




.53 


.AAr 


.57 


yoLun. coEF. 


.002LA^ 


.002G 


.OOZi- 


- 0 033 


VFlocity 


Z.0.O 


ZO.o 


Zo.o 


20.5' 


FUEL VvO", 


I.9T8.7 


2,378.8 






H/1CH. WT- 


8H3.G 


ssr.i 




?7o.8 


OUTFIT WT. 




3,‘ioS’, / 




3.V/6.^ 


STEEL VVr 




ni. i 




6S7S-0 


Mist, WT, 


3oo . 


300. 




300. o 


Ballast wr. 


ze.i 1 i.H 




2, /4V,*7 


Pax load WT 


3’,T65,3 


?, (=s9.o 




/Z, 22^.0 


• &M 


IHS 








Klo. CO/rrA^^/£f?S 


Zio 


58sr 




IHO 


LEV. roT^i-D^ 


0 


\ 




Z 


LEV. t^cOC 


3 


£T 




5" 


No. Shafts 


Z 




\ 


1 


No. FajG-ine5 


L 




z * 


Z' 


SH? REd. 


'3, 


/L,?/7 




21,000. 


COST 1 


.io.o 


2 2.0L 




Zh.S\ 


AOC 




tr. Co 




G.S\ 


RFR 










CRF 










FR 


ICO. 


/oo. 


loo. 


loo. 


WPV 


(=) 


Z/.^tT 


© 


27. S’ 


?RofeMFO 


\k)7=>u7 


d:>c4rPtJtr 


(>UPU7- 


<outPut 



* Fo/^ce't) cxtRa For Reci/^^fciTV. 

|^T>iCATH^ WFG‘A'nv/f^ \/ALVB 

Table 7 (cont.) 
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counteract this tendency. Suggestions for future modifica- 
tions are covered in the concluding chapter, 

9.4 STOiT Trajectories 

The following tables show the trajectories followed 
during the search by the SUI.IT program. The starting point 
and program optimum alternatives are given in terms of the 
decision variables. The intermediate solutions are shown by 
the incremental differences resulting in the changes between 
successive steps. The tables also indicate which constraints 
are constraining the solutions. These intermediate solutions 
represent the solution for each of the subproblems associated 
with the successively decreasing value of the program para- 
meter 'r'. 

Most of the improvements in the objective seem to occur 
in the f irst f ev/ iterations. This indicates that the process 
could be streamlined by imposing looser convergence criteria. 
This would permit the program to terminate its search earlier. 
It v/ould not be necessary for the penalty terms introduced 
in the modified problem to he diminished by an even i.urther 
reduction in the parameter 'r'. The execution times for the 
optimization procedure are given for each of the problems. 

The comparison between results is not possible here. The 
differences in the objective function result from various 

freight rates used. The final table of comparisons account 
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SOLUTION TRAJECTORY 



DESIGN 

VARIABLES 


STARTING 

POINT 


CHANGI 

SUBI 

1 


E IN ' 
PROBLEI 
2* 


\/ALUE 
M NO. 
3 


IN 

4 


FINAL 

POSITION 


Displacement 


So, fas'?. 


f 0,‘>9o. 








39, 


LENGTH 


ns^o. 


-11 








^2.3. 


BEAM 


93. 


4X. 


t .1 






9^.1 


DRAFT 


2’?.3 




+ .9 








DEPTH 


55“. 


•+3. 


4.n 






^8n 


PRISM. CoEF. 


.tr 




— 






.492 


velocity 


-2?.^ 


-Z.’B 


-.1 






72.8 


CONSTRAINT t.mc 
VALUES (Sfc) 


I.OI o 


12.390 


H.OTO 






14 . oS 0 


freeboard 














STRENGTH ^ Vd- iS”) 














Max. length (9po) 














Max. Beam (no.) 














MAX. draft (37) 














Max. %: (i.a) 














Min. Vjl (•?) 














MAX. CP (.70) 




TiG4a 


Tl^yNT 






7 ((tHT 


MlN. CP (.98) 














MAX . % (^ vs") 














M/N. Vr 














MAX- C.V i.oob) 














MIN. Cv (.001) 














CB ~ .9£5 CP 




TIG-hT 








7'iG'HT 


OBTECTIVS 4 
( K/PV X lO***) 


3F.8 8 


+ 89.0 


+ 1. 






129.88 



+ SouirrioN At C^) Iw’tiCA.T^T w% CH4A»0f 

OF 5W?^P>?oBLfAl *2., 



*■ FR-'^ZOo./y^ 



Table 8 
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SOLUTION TRAJECTORY /^.ternative^Z® 



DESIGN 

VARIABLES 


STARTING 
POINT • 


CHANGI 

SUBI 

1 


E IN ' 
PROBLEI 
2 


\/ALUE 
V NO, 
3 


IN 

4 


FINAL 

POSITION 


Displacement 




+ 5 4(,0 


— 


■» (0. 


- 


45,310 


LENGTH 


“123.0 


^ 38 , 


- 


- 




7bi, 


BEAM 


*?‘t.| 






-- 




?e.5 


DRAFT 


32.2 


+ I.ZC 


— 


— 


— 


33. H 


DEPTH 


5-«.7 


4 .g4 


— 


— 


— 


59.r 


PRISM. CoEF. 


.(.n 


— 


— 


— 


— 


.(,92 


velocity 


22.8 


+ .S 


— 


— 


— 


Zl.l 


CONSTRAINT tim? 
VALUES CsfO 


.960 


A.A60 


S.^o 


6.46p 


0.270 


8,270 


FREEBOARD 














STRENGTH ^ Vd- 1^) 














Max. length (9oo) 














Max. Beam (mo.) 














Max. draft (37j 














Max. %: (i.e) 














Min. V/l (•?) 














M/IK. OP (.70) 


TiCtHT 


TlCH^T 


T'\6,wt 


TIG-h-T 


TI^HT 


T(CtVv7 


MIN. CP (.G0) 














MAX . ^/t (3-75) 














M/N. (^-^) 














MAX. <^v (,o(?fc) 














MlN. Cv (.001) 














CB ~ .925 CP 


T(6rV+^ 


W-T 


T< 0 


T(GrU7 


TrCH-T 




oBTEcnvS * 
( K/PV X 10'*“) 


/ 

12t.2 




-to. 


to. 


t O. 


137. fe 



PA '•^200./t-c.w 

C-\ WUKATpk A>« CH<^wGtr Table 9 
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SOLUTION TRAJECTORY 



alternative 




DESIGN 

VARIABLES 


STARTING 

POINT 


CHANG 

SUB 

1 


E IN 
PROBLE 
2 


VALUE 
M NO 
3 


IN 

4 


FINAL 

POSITION 


Displacement 


1 s;ooo. 


>•7,140, 


-1840. 

1 


*2 ft 80 . 


-9. 


/2,7^|, 


LENGTH 


boo. 


o. 


+ 1, 




— 




BEAM 


So, 


4- o 


- 1.9 




— 


7g.Z 


DRAFT 


15. 


■*s.z 


-1.8 


+ O 


— 


23.<< 


DEPTH 


"iS. 


4 O 


+ 0 


o 


— 


^5-. 


PRISM. CoEF. 







- 


~ 


— 


.to 


velocity 


Zo.o 


- 


— 


- 


— 


-2<?.o 


CONSTRAINT t/mc 
VALUES (s-cc) 


.960 


T.b^O 


10 .do 


12.500 


/3.370 


13.370 


FREEBOARD 










TvOmT 


TltrH » 


strength ( %- 










TIGWT 


TUj-hT 


Max. length (9oo) 














Max. Beam (ho.) 














Max. draft (37j 














Max. %: (i-z) 














MiN. ('5) 














MAX. OP (.70) 














MlN- CP (.^0) 














MAX. ®/r (3-v^) 














M/N. Vr U.is) 














MAX. (^v (.oob) 














MlA/. cv (.001) 














CB ~ .9^5CP 














OBTECTIVE » 
{ KJPV X 10**“) 


l.Sb 


+ E7.I 


•V l.t 


+ 5-.Z 


4 0. 


■J1.M 



# 



FR^ 



200. /ton 



(-) wDKArrs MO Table 10 
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SOLUTION TRAJECTORY alternative©-® 



DESIGN 

VARIABLES 


STARTING 
POINT • 


CHANGI 

SUB 

1 


E IN ' 
PROBLEI 
2 ^ 


^ALUE 
M NO, 
3 


IN 

4 


FINAL 

POSITION 


Displacement 


1 E, ooo. 


+ 13,38V. 


— 






20,38A, 


LENGTH 


boo. 


■*'10. 


— 






670, 


BEAM 


?0. 




— 








DRAFT 


2.F. 


4 7.b 


— 






32.4 


DEPTH 


^5" 




— 






^9,7 


PRISM. CoEF. 


.4 


4 .08 


— 






.48 


velocity 


2o.o 


+ .5 


— 








CONSTRAINT 

VALUES 














free-board 














STRENGTH 














Max. LENGTH (900) 














MAX. Beam (ho.) 














Max. draft (3t) 














Max. Vrc (i£) 














Min. (•?) 














M/liL. C-P (T") 














MIM. CP (.^8) 














MAX. % (^VT) 














M/N. ®/r 














MAX- (,o(?fe) 














HIN. Cv (. 001 ) 














CB - CP 






yiir^T 


« 




T/(r 141 


OBTECTIVE , 
^ NPV X lo***/ 


4.48 




t ,0 






54.32 



0. <L^) ^‘^tXCArcS AJO CMAV6-V 



+ SocutioaU CO/ove^Wl AX E^X> of 

sub?/ro'bi.fh *~z. Table 11 
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for these differences and presents the results based on a 
freight rate of $100/ton, 

9.5 Ranking the Alternatives 

After conducting the individual optimization problems 
with the various starting points, a comparison of results was 
desired so that the best design could be identified. The 
various alternatives v/ere compared using three measures of the 
economic worth of the designs. The following table gives the 
Net Present Value of the design based upon $100/ton received 
in revenues. Also the Capital Recovery Factor based upon 
the total investment and the Required Freight Rate are listed. 
The RFR depicts the rate that must be charged by the opera- 
tion of the ship so that the revenues would just cover his 
costs which include depreciation charges. 

The matrix shows that no single ship is preferred using 
all these criteria. It also shows the improvement obtained 
in the measure of effectiveness using the SUMT Program. Al- 
ternatives 3 and 4 switch position in relative ranking as the 
objective is switched from NPV to CRF or RFR. The solution 
to the scenario with the given assumptions would be Alterna- 
tive 4, Other scenarios may result in one of the other de- 
signs being a solution. Each of the designs carry a different 
payload. The differences in worth per ton transported is not 
.so great which indicates that two smaller ships may be 
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A 



I 




RANKING OF ALTE-RNATIVES 
BY VARIOUS CRITERIA 



alternative 


NPV 


ERF ^ 


KFR ^ 


A- 30,659. 

Viv j- Shaft 


10.48 [s] 


,13 d] 


82.37 5 


^ ^. 30 , (.59. 

\.£y Z Shafts 


(-) □ 


(-) □ 


loo-tr B 


A= 39,64-1. 1 
VRY 2 shafts I 


ks.os H] 


,7f □ 


35.35 ■ [D 


A ^ ^s;5io. i 

VlY A SHAFTS ! 


[5i.?o DD 


.57 d] 


43.70 B 


A - 15 ^ 00 . 

N ' 2 S HA FT5 


(-) □ 


(-) □ 


1E3.T9 0 


V_y 2. SHAFTS 


K-TS- [T] 


.37 B 


47-9A B 


4 = /5',ooo, 
V-V X Sh/;ft 


(-) □ 


(-1 □ 


- □ 


A =^,384. 
N — A 2 Shaft 


27.10 d] 


.^8 II] 


42.80 B 



E IMTJlCATFi ■REtAT.o'f RAKM-Mfr 

(/JblcATeS «f6/tmuF VALUe 



Table 12 
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preferred over one larger ship. The increased flexibility 
available to the owner may make a number of small ships an 
attractive proposition, 

9.6 Discontinuities 

During the design of the various alternatives, discon- 
tinuities were observed in the objective surface. These were 
the direct result of the discrete characteristics of the en- 
gine and installed cargo. The following graphs show the 
effects of the discrete factors in the design and additionally 
shov/ the performance of the search routine in the neighborhood 
of the discontinuity. In order to demonstrate these factors, 
the ship hull design associated with design Alternative 4- 
was held constant as the design speed v/as varied from nine- 
teen to twenty-seven knots. This provided a profile of the 
objective surface along a single dimension. A similar graph 
is also developed for Alternative Number 3. 

The effects of the discrete power plants are evident as 
the speed is increased. The graphs note the speeds at which 
additional engines are required. The cost function reflects 
the additional costs involved in the transition. In the graph 
for Alternative 4, the quantum drop in the objective function 
just before twenty-four knots is due to the cargo capacity 
lost when the engine box was enlarged. The same integer num- 
ber of containers could no longer fit in the length available 
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OBJECTIVE SURFACE PROFILE 
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and this resulted in a loss of about fifty containers. As 
stated earlier, the length of the engine box was modeled 
solely as a function of the required shaft horsepower. If 
the engine box were modeled more accurately, the loss of cargo 
may have occured at the same time that additional engines 
were added. The graph of Alternative Number 3 shows a simi- 
lar drop in the value of the objective function. This is 
again associated with the engine box size. At the bottom of 
each graph, the tight constraints are listed for the various 
design speeds, 

9.7 Sensitivity 

The final segment of the analysis is concerned with 
determining the sensitivity of the design solutions and the 
associated objective values. Only a limited effort was pos- 
sible in this area, but this section should help to re-enforce 
confidence in the results of the previous sections. The re- 
sults of this step will point out the critical assumptions 
made in the design model and v/here additional effort may 
prove rev/arding. 

In order to obtain a sensitivity analysis in a general 
search routine, it is required that perturbations be made in 
the values of the decision variables and the input parameters. 
Here we will study the sensitivity of the design Alternative 
4 due to variations in discount factor, propulsion coefficient, 

107 



midships coefficient, endurance , freight rate and the number 
of extra engines that may be required for increased relia- 
bility. 

The procedure was to start with the best alternative 
available, and determine the affect of changing one factor 
at a time. The clianges on the worth of Alternative 4 are 
noted for each perturbation. In Column B of the table, the 
NFV objective resulting from the change is recorded. This 
design was used as a starting point for the SUI^T Program. 
Column C indicates the v/orth of the design resulting from 
the optimization. A qualitative indication of the overall 
changes are also given. 

The most critical factor in the design v/as the assumed 
value of the discount factor. A variation in this input 
parameter caused the greatest change in the design. There 
are several other elements that deserve explanation. 

The sensitivity of the design to the assumption of P,C. may 
be more important than indicated by a quick look at the fig- 
ures. The reduction in the P.C. resulted in a design with a 
higher block coefficient. This would cause an even further 
decrease in the P.C. for the design. The other interesting 
effects are noted on the table. The design incorporates 
a severe cost on the design for ballast, V/hen the ballast 
is replaced by fuel, there is a sufficient increase in the 
NPV of the design. This appears twice in the analysis. 
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CHAPTER X 



CONCLUSIONS AND RECOI.H>^.ENDATIONS 

10.1 Genera-l 

As the first iteration of the design procedure is cul- 
minated for the sample problem, it is important to put the 
discussion in context with the overall design philosophy. 

The preliminary design was pursued in this problem v/ith a 
broader interpretation of the owner's requirements as was 
made possible by the quantification of the owner's prefer- 
ences. It was found that this approach resulted in a 
streamlining of the design step which permitted a larger 
range of alternatives to be evaluated. The results from the 
sample problem, would be presented. In the proposed design 
process in figure the next step would involve an evalua- 
tion of results of each system design by the owner. This 
would be follov/ed by a possible decision to iterate in the 
design of the ship with modified criteria. 

The ourput of the design step should provide useful in- 
formation of the ov/ner so that future iterations will con- 
verge faster to the design selected. There will be many 
occasions where the owner will want to have results of a 
design based on the traditional requirements of speed and 
payload capacity. This approach used should prove adapt- 
able. The information on sensitivity would be obtained. 
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10.2 Optimization in the Design Process 

We have seen that the computational efficiency of the 
computer has permitted extensive parametric study of the ship 
design. The interaction between the problem formulation and 
the solution technique has been demonstrated. The application 
of the Sequential Unconstrained Minimization Technique re- 
sulted in a straightforward design sequence that could be 
used to evaluate many design alternatives. The ship design 
model was developed to provide the information necessary 
to identify the objective surface for an optimization proce- 
dure, In the sample problem, the surface was generated using 
the net present value for measuring the cash flows expected 
over the ship's lifetime. 

The thesis has shown that partial optimization of a 
ship design is possible if the problem can be structured. It 
was also observed that interaction of the designer with the 
program is needed to achieve the maximum benefit from the 
programs , 

10.3 Experience with SUi'.'T 

The experience with the SUIv!T program generated several 
observations releveant to the design process. The first is 
that a programmed search v/hich requires derivatives of the 
surface, could be employed using linear interpolation of the 
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derivatives and remain relatively efficient. The SUTflT Pro- 
gram used such a linear interpolation scheme to determine 
the partial derivatives of the objective surface. This 
efficiency was improved by deriving the design constraints 
so that interpolation of their derivatives would not be 
necessary, Por those constraints which were not expressed 
in closed form, they v/ere handled internal to the design 
model by insuring that the ship designs would satisfy these 
constraints. The non-linear stability constraint was handled 
this way because the change in stability could not be deter- 
mined from the decision variables along, but must include 
the effect of the container loading. The design model was 
constructed to insure that all alternatives would meet the 
stability constraint by controling the loading of containers. 

A second observation dealing with data restrictions 
should be imposed external to the design model, V/hen they 
were introduced internally, it was difficult to assign 
appropriate penalties within the design model to keep the 
design feasible. No matter what reasonable level of pen- 
alty v/as assigned, the program v/ould find alternatives that 
did not satisfy these constraints. 

The imposition of the equality constraints was not 
found to be practical in this design application because 
none of the intermediate designs were feasible. It was only 
after the parameter "r" was sufficiently small that the 
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equality constraints imposed sufficient penalty to the modi- 
fied objective function. It was preferred that the equality 
be replaced by only a single one v/ay constraint. If two 
inequalities were imposed, the associated penalty in the ob- 
jective function would be of infinite value and as such v/ould 
obscure the real driving factors in the problem. This 
approach was used for the equality constraint of the block 
coefficient in the sample problem. The nature of the prob- 
lem ca\ised the program to design ship alternatives with as 
large a block coefficient as possible. The equality was 
replaced by an upper bound. 

There are two other observations on the application of 
the program v/hich should benefit future users. The first 
deals v/ith starting points for the program and the second 
with the selection of the value of the parameter "r". 17ith 
the regard to the first, in the present form of the SUI'T 
Program, it is recommended to start the search with a 
feasible design. The program is capable of finding its own 
feasible starting solution, however, the barrier penalty 
technique used in the optimization program requires a sepa- 
rate method for finding a feasible starting solution and 
attempts to employ its use has resulted in imreasonable 
alternatives, A feasible starting solution v/ould prevent 
the excessive costs experienced when this feature v/as tried. 

Finally, it is recommended that the selected value for 

"r" and for the associated reduction ratio which determines 
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the successive values of "r" be determined as follows. Look- 
ing first at the modified objective function, it is impor- 
tant to have the values of the constraint and objective 
functions of the same magnitude during the first iteration. 
This will insure that the program encounters both the affects 
of the objective surface and the constraint penalties in the 
initial phases of its search. If the value of the constraint 
were larger than that of the objective function, the first 
several iterations would produce solutions that did not ob- 
serve the contour of the objective surface. If on the other 
hand, the objective value were larger, the program may sat- 
isfy the given convergence criteria prior to reaching a 
local optimujTi on the surface. 

In the sample problem this sizing was accomplished in 
tv/o ways. First, the objective value was expressed in 
millions of dollars. This resulted in the objective value 
being of the same order of magnitude as most of the penalty 
values that were possible from the constraints. For this 
reason, the value of "r" was set equal to one. Other com- 
binations have v/orked as well. The selected value of the 
reduction ratio depends on the convergence of the sequence 
of subprogram solutions. Most of the improvements in the 
sample problem optimizations v/ere made in the first few 
iterations. The remaining iterations were to decrease the 
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penalty values so that the convergence criteria would be 
satisfied, I would recommend using a higher ratio than 
twenty which v/as used in the sample problem. 

10,4 Areas for Additional Study 

Future applications of the design program could pro- 
vide even more realistic results if the following improve- 
ments were incorporated. They were not included in this 
prograun due to the limited resources available. The follow- 
ing comments cover most of the areas that could use addi- 
tional effort, 

Fro-pulsive Coefficient 

The sensitivity analysis showed that the solution is 
sensitive to variations in the propulsive coefficient. As 
the design alternatives develop blunt sterns, the efficiency 
of the propulsor will decrease. The ship propeller inter- 
action should be included. The determination of the actual 
propeller would also prove helpful. One approach to deter- 
mining propeller fit is covered in Reference 45. 

Other Constraints 

Additional constraints coimnonly handled, such as the 
requirements on the period of roll should be included. This 
would require only a slight change in the program. Also 
the volume requirements v/ere used only to determine feasible 
container configurations. The addition to the model of 
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restrictions due to fuel and ballast volumes should be de- 
veloped internal to the model, The volume requirements of 
the engineering space should be more carefully defined, es- 
pecially since the length of this compartment is used to cal- 
culate the length available for containers. 

Velocity in a Seav/ay 

The additional resistance experienced in the seaway 
causes a loss of speed. The program uses a simplified model 
to predict this loss. A subprogram would prove useful in 
determining the value more accurately. 

Interaction v/ith Terminal 

The program presently uses a very simplified inter- 
action measure. The model could be expanded for a specific 
situation to determine the actual turnaround times and 
charges. This v/ould require information on the container 
handling capabilities, both in terms of speed of operation 
but also the inventory capacity. 

Backhaul 

Any actual study should investigate backhaul opportu- 
nities and incorporate their effects. 

Maintenance 

Maintenance was assumed to impose an annual dollar 
cost. If an overliaul policy is knovm, the cash flows should 
be adjusted. The Net Present Value objective function is 
the only measure v/hich could handle the non-uniform cash flow. 
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Seakeeping 

An effort should be made to determine the seakeeping 
qualities of the design. At the least, limits of perform- 
ance could be included in the constraints. One such limit 
would be the deck wetness criteria of Reference 39. 

Pov/ering Prediction 

This program used the Taylor Standard Series to predict 
ship resistance. Other predictors such as the Series 6o 
could be easily employed. Of primary concern is the power- 
ing predictions for designs with large block coefficients. 
Another alteration would provide calculation of resistance 
for the different ballast conditions. 

10.5 Applications 

The design model was developed for the preliminary de- 
sign of a containership. There are other ship types which 
could be studies with the existing program. For example, a 
LASH or Barge Ship could be considered, also bulk carriers or 
oil carriers could be designed. This would require an ad- 
justraent in the input parameters describing the size and 
weights of the containers. The" barge is a large container 
and the analogy is obvious. The bulk carrier on the other 
hand is loaded in a continuous way. By loading small 
blocks of oil, the program could approximate the cargo 
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carried. 

Another application of the program would call for its 
use to test the effects of governmental policies on the 
nation’s merchant fleet composition. The effect of subsi- 
dies has the effect of changing the economic preferences 
of ship owner’s for different ships. The model could show 
the trends in the fleet population that were due to these 
policies and regulations. Of more immediate concern is the 
effect of oil prices on the design. 

The last suggested application involves using the 
programs as a teaching aid in a graduate course in ship 
design. The following characteristics make this applica- 
tion attractive. First, the student is required to provide 
a feasible starting point for the design. This could be the 
result of a hand calculation and as such provides an oppor- 
tunity to understand the design methods. Secondly, the 
design program could then be used to generate the data to 
check the student's design accuracy for his selection of deci- 
sion variables. Finally, the optimization program would pro- 
vide other design alternatives so that the student could be 
exposed to a higher level in the design process. As such, 
this would provide an interactive tool to permit a more 
complete design effort in the more advanced design courses. 
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Appendix B 
EQUATION DB^INITIQN 



Contents 



Equation definition for design model. 

Graph of speed reduction factor for speed made 
good during transit. 

Curve of machinery weight as a function of shaft 
horsepov/er . 

Curve showing model for determining SPC for gas 
turbine pov/er plant as a function of shaft horse- 
power. 

SFC plot for various gas turbines. 
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Empiricol unci Other Equotiont Used for Corgo 
Ship Modol 



1‘lif |*n*M‘iiU*il in lilts ;i|i{U‘:uli\ ;irt* 

tuken from ri-tVri iiiv 1 1 j iiihI I'onii the l>;ivis I'nr the 
c'uiii]>nt:itioiis tli;il ;irc eiirru'il uiit in >U*|i •! o( 
cncli simpliu;; eyrie. These e<|n;ili«ois cmiM 
uiKlouhtully he presruleil in iliOVrent form with 
im|>nive<l ;ieeiiniev. The ei|ii;itious are lisleil in 
seven major r.roii|»s: ( I ) \\'ei;;lils. 11*. (21 Voinmes, 
K, (3) ITeehoanI, /*', (I) Stability, 6*. (a) Cost. C. 
(0) Power, f*, mill t“) Misri'llaiieoiis. \I. 

I Wa'iififx (all weii;lits are .I'ivt ii in Iimij; tons). 

(a) Approximatiim t<i onttit wciKHt 

Kt O.OSO 



(b) Appruxiniutioii to steel weight 

07 ['i? BUI 

(f) . Appro\iuiation to wet machinery wciijlit 



W, - 7.1S SHP* *“ (3»0 



(if) Weight c,f lucl oil required tu sail tne 
specified distance. £ plus 10 percent alloivance. at 
a given speed. (1"). 



W, 



(l.lll£) XSH P X FK 
2240 X V 



o.4!)i X nr* 



(c) Approximate fuel rate 

FK - 0.5 SHPASlIl* - .S.V,) (5 ID 



(J) Miscelluiieousileadwciglit 

W, -300 (OIF) 

it) Actual payload weight 

\\\ - IP. - II'„ - IP, - IP, (7110 

where A » .Vi (see Table I). 

2 VoUtmes (all volumes given in cubic feet). 

(a) Approximate gross hale eubie of the ship. 
This value applies to all dry-eargo .ships exeept 
where e-xco'^ivc sheer is use<l. This expression 
includes machinery space volume ami excimics 
thedouble-boltoui, settler, and peak tanks 



GlIC « O.S75(/u XBXDX C\I/v% (IT) 



(h) .\ppro\inia(r I'm I oil rapavity in llu tloiihli- 
Ihitliiin. rhi> Kpialioii i^ used in roujiimliiat 
with e(plati*iit I *) loiU ti rinitte u lirtlii r adripiate 
vohitiu* exists to arroiiiiiiodate the fuel oil coin* 
pitted by ei|ii:iti<m t ( IP). 

(VoI)„ - ((A\ X L) X/iX (A, X D) 

X (O.G'J C/,)I 

« ii.7N.SA'. XGHC (2P> 

where 



A', es p. n 

C'k * hhx“k ciK'llieieiit at LWL; factor b‘.) |kT“ 
cent is u eorreetatii for (o) sirnctttre in 
inner boUotii. and (A) correction to obiaiit 
Cn at the \VU height eiptal (o tank top 
height. Puel-oil stowage factor i> ^17.2 
cu ft/ton 

(<;) Approximate furl-oil seUler-Uttk capacity 
(Vol)/r - oOlH) cn ft d.*)0 tons) (:iP) 

(d) Approximate nmcliiticry space vedume 

(VoI)„ - -I7.0.T> -f 7.112 (SUP) (-IP) 

(e) TIic actual payload volume is 

Vol, - GBC - VoUh,..e, - (37.2 IP, 

-(VoU-f VoV)) (.jP) 

\j / A ilC 

SF - Vol,/IP, (GP) 

3 Freeboard (all freeboard dimensions are in 
inches). 

(a) Available freeboard. Equation assumes 
3*in. margin line and that the uppermost con- 
tinuous deck is the freeUiard deck. Equation w ill 
differ for a shelter-deck ship 

F. - 12 (/) - (0.25 -f T)] (IF) 

(b) Mtitimuiit |)crntissible freeboard from 
USCG Loud Line Uegnlutioiis. 

For: A ^ -HH)ft: 

£, a 4.21+ 3.50 -|;; + 3.71 (2/0 



For: -100 ft A ^ 750 ft 



Fr ^ -77.07 -4- -12.5S 



lUU 



FoK L2. 




(• 1/0 
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4 ShihUlfx 

(«i) Apiniisiiniiiiun P‘* tin- tr,iii>vus<- iinrti;i 
ill’ llu* lu.hl \v;tt«r |il.iili*. <f 

tr.ms. w:iti‘r|il.iiic uu'Uiu/Ui' 

a » (Min;»7 C‘, - 0 .MIJ 2 (15) 

(A) A|>|>r(t\iiniitMiit l«i thi* f\(! of stivl 
K(*t * K% i) 

K% - iUn (25) 

(c) ApproNiimitiiHi tu tho KG of oiiUlt 

KC. - A', n 

A% <• {AM) (35) 



uliifh ;iiv rt l:iU<l to <l«ill:ir\ ity a Victor i>f 
lirlmcii I ;hhI ~i\ '' 

(li) OiiUit co-sts 

Knr; 0 i ir, S I Mki pmis 

(; » |inM) - UJM3 n\ H (0.M2 X m*-*)ir«- 

- (<1.1323 X in ")IIVl 
Tor: Him < U‘, < 2i^m phin 

C. - 12430 - I.02S \\\ 4* (0.722 X IO->) n'.» 

- (0.001 X 10 *) nv\ n\ (20 

For: ir, £ 2000 tons 



(d) Approxiuintimi to the KG nf psiyloucl 

KG, - A*7 D 

A% • 0.(VJ (45) 

(e) Approximation to the KG of ini*;ccUaiicous 

deadweight 

AC7, - A, D 

Ai « 1.00 (35) 

(/) Approximation to KG of fuel oil in settler 

tanks 



KGf, - A*/) + 4.S0 

A, - 0.11 (65) 

(f) Approximation to KG of fuel oil in double 
bottoiO Prr f.rr**:*'* rW.'p- (iV^) 

(tank top height) 

KG„ - 0.07 K,l) (75) 

(A) Approximation to machinery AG. with 
boilers full mid for convent jonal arrangement for 
steam-turbine plant 

AG, - 0.53 D (S5) 



(i) The available CM in the full load condi- 
tion uncorrected for free >urfucc is approximated 
by 



-((ir. KG, + \l'.KG. + ir.AC, 

+ (II', - W’ft)AG,. (!15) 



+ H'„ KG„ + ir, KG, + l^, KG,]/S 



C, - 0!).S.S II'. (.10 

(4) ApproNimation to steel cost 

c. - [2I.V4 - •->l.•^‘=(^;;-;) 

+ 2."ci(^,y-o.iu9(-;^^)’] It'. (40 

(c) Approximation to machinery cost. (SHP 
is normal shaft horsepower.) 

• For; SHP ^ 13.IKI0 



C, - 



137.7 — 



SHP 



75..32 4- (SHP) 



SHP 



For: SHP 2; 13.000 

C, - 



r SHP 1 

3.249 (S14P) - 173.05 ’ 

L 100 



SHP {DO 



(d) Annual fuel cost 



Cf - A,oAMIVI7£ 

Aio ■» 5205 for this study (7G) 

Ai — 3.09 cost pls./ton t 



An explanation of A'lo is given in section S of this 
Appendix. 

(e) Total annual cost. Cost 
The yearly fuel cost (c<iuation (70 1 is included in 
this cost which is computed as follows: 



Aj » 0.54 

(Jj The rccpiircd GM is 

G.IA i 0.05 li (105) 

(This rctpiired (V\I is also for the /full-load condi- 
tion, uncorrccicil for free hurlnce. j 
5 Cosfs (all costs are given in terms of cost 



Cost * 0.UI550 X (C. + C, 4* CJ 4* Cf (SC') 

where the fraction O.IMmO’ reduces the total initial 
building cost to a yearly cost using the assunip- 
tions: 



’*Thc Urm C(»l points is uvd in (1| to avoid direct 
assnciatiitn with doKurs. 
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REDUCTION FACTOR FOR SPEED MADE GOOD 
DURING TRANSIT 




Figure 19 
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Figure 20 
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SHAFT HORSEPOWER 
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Figure 21 

135 



SHAFT HORSEPOWER (X 1000 ) 



EXPLANATION 0^ Sl^C MODEL 



The graph in Figure 22 shows the minimura 
SFG as engines are added to the power plant to 
achieve a given installed shaft horsepower. 

The initial segment shows the SFG curve for a 
single engine. For the multiple engine install- 
ations, the same minimum SFG is achieved at in- 
stalled shaft horsepov/ers of multiples of 18,700 
SHP. For the intermediate powers, all engines 
are assumed to operate at equal levels. 
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Figure 22 
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SHAFT HORSEPOWER ( X 1000 ) 



°/o SFC 




SFC CURVES FOR MARINE GAS TURBINES 



Figure 23 
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N.N.S. DESIGN 



PROGRAM INPUT PARAMETERS 



PARAMETER 


SYMBOL 


VALUE 


DPSl&M ?/^/^/\MGTfcRS 






PROPU LSit^e: COEFFiClGAJT 


PC 


.Ti (aW>P0M6-6 1-03) 


EfJX)UHAhJ Ct 


E 


lX.Soo. 


D^S7/l^Jcf BertJeE'^ PoRtb 


DBP 


SOOO 




A-VA 1 L 


■9 8 


OP 


nsha fT 


1 


■Ml^^. NO. OF 


NGAJGti 


(-:> 


wo. OF SP/'RF 


AyexTRA 


o 


HA)<. pAVcoAb WT. 


WW 


Xo^ooo 


MW. P/^VLOAb WT. 


V/PMIA^ 


O 


Max. COa/tai/^SA S7/}CK(A/6- 


KJV 


'T 


MAX. STACKIA/O ABotA' t)K. 


MD 


1 


AVB/^A0>F COa/TAWcR vJT. 


WC 


H.8G 


BwecnuS <!oajT. lsaigth 


ELENCrT 


33.0 


ETTPOTIIXE COMT. WiTrrH 


£ WfDTH 


8.7S 


BrFeaioe t)f'FTtf 


Et>FPTH 


8.02 


T)0((Blf Bottom TH»ckWEis 


t)ou B L-P 




BALLAbT wet&t-ir kG. 


K6-^ 


:2.o 


Max. 5HP Fo/X Si'yj6<.e BAJG. 


Sm?(A 


N.A. 


SSA 5Pt=fT) FActo/P 


BftA 1. 


■ OO'TS" 


5 PC. FACTOR 


Beta 2, 


•34F 


Ml SC. WF/Fivr 


wxx 


O 


Mist. w/griG-^rp K&. 


XK6rX 


O 








C>E5I6a; vVE-((B4Pa;F factors 




- 


Factor, for ouTF(r irt. 


. WWo 


1 . 


FAC.TOR f<3R stpfc wr. 


WWS 


J , 


• Factor for ma<h<a/£A‘;' va 


S</\a/k1 


1 . 



^ VAL^i? ^TiFS TUAf^tAJir Pf^oPUL Si 0/0 , 
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N N.S. 



PROGRAM INPUT PARAMETERS 



PARAMETER. 


SYMBOL 


VALUE 


NAVICATIOM 






MAX. LEM&TH 


LEa/GTM 


9oo 


MAX. BEAM 


13fcAM M 


llo 


MAX. tXAFT 


DRfl FTM 










Ecowom'C. PAMtArrefis 






PRFiCrHT rate 


FR 


so. 


<LALE F££ 


PCF 


looo. 


PAiLY 'Pofi.J 


T>P F 


looo. 


PRlCt Of- Fa^L 


PRICE 


lEs 


AWVUAc iMSURA^Ce P^AIC 


AXR 


. 1 


A5SCT El Ft 


NN 


IF'. 


SALVAGh VALUE 


SV 


O 


Discount factor 


DF 


.2.0 


’•SuBS'T>V I^ATF 


SUS IDR 


.SE 


PBiRQjF^ Fl>^ANCT-'t> 


PF 


.73' 


\HTBREhT RATP 


RATE 


JO 


CoiCfP. Tax RATF 


CTR 


.^€ 


iNVFSTMFA/r tax CAebIT 


AfTC , 


.07 








C’OST W<TlGHTl^G factors 






FAcTO/^ Fo/l OJTFtl COST 


CCo 


1. 


Bk-to/^ fo^ 5T&ec Cost 


CCS 


1. 


FA'CTOA. FOA. MAQ4(A;e<Y c.. 


CC-M 


t 








PRoGMM PARAMETERS 






OBjFcnoc Sccectioa/ 


l<pP>l 


I (ajTu^ 




DcLtA ex') 


Ate Co ) 



Table 15 (cont.) 
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B.I.W. DESIGN 



PROGRAM INPUT PARAMETERS 



PARAMETER 


SYMBOL 


VALUE 


DP516M P>A»^/qMn-eR5 






pRoPu LS(t'cr cogffic.(ga;t 


PC 


•To (ATT>eATM6e l.os) 


EKiDWRA hj tt 


£ 


12, 500, 


D(STAf<)ce Betw<=ea^ Ports 


DBP 


5^000 


Si4jT> >AVA'L/1BlL<rV 


AVAIL 


.9 0 


OF shafts 


nsha ft 


1 


hlF/. NO. OF EAj(rt/UFb* 


HBfoeiJ 


(- > 


fJO- OF S'PARF eW6-(>OFS 


A/GXTRA 


O 


Ha/-. pA'iiOAb w/T- 


WW 


Xo,ooo 


hiiAJ. "Payload wt. 


WPMIAJ 


0 


Max. <IOa/7/1IA/Ga! S7/4CK(-^6- 


Mt) 


3 


MAX. stackia;o ABo</t dk. 


Mt) 


6 


AVFi^AOE' COaH^/VcR. h/T. 


we 


1 \.T 


BHC-CTIUS CO/OT. L£A)0rTli 


ELEN&T 


23.0 


eW€<FX\K/0- COKiT^ VVI177T4 


£wn>TH 




BrFBCnoE Coa/T. X>B'PTlf 


EbFPTR 


S.OZ 


DouBlf Bottom TH«:KA;ois 


t)ouB,cp 




B/lLLAsr WeiANT NO. 




2. .S' 


Max. 5H-P FoK S(Aj(u.e fajc,. 


S|4?M 


M.A . 


SF/A SPffT) FActo/P 


Bft/? 1 


.OOlF 


Spa pac-toa 


Beta 2, 


.365' 


Ml sc. WF-ldrWT 


wxx 


o 


Fuse. uyFfG-HT KG-. 


XK6-X 


<9 








DfsiIsaT u/fc(644-nv6- FAOTO/^S 






Fa CTO A. FtP>^ OuTFrr 


\</Wo 


(. 


factor for stffc wr. 


wws 


1. 


• Factor for MAceH/OeAYw 


w/ki 


1 . 



^ VACM^ l'Ot>lC/^TBS ST<r/^ TU/<0 ^a^C- P^OPUL^tO^ , 

l42 Table l6 



B, 1. W. («o*n^ 



PROGRAM INPUT PARAMETERS 



PARAMETER. 


SYMBOL 


VALUE 


NAV|C/\TI0N CoUil 






MAX. LEM&TH 


LEA/GrTM 


90O . FT. 


MAX. BEAM 


■beamm 


l\o. ft. 


MAX. D«AFT 


DR A FTM 


5D. FT. 








ECOWOMiC- P/>7V\MET£-R5 






FRPlCr(4T RATF 


FR 


5^0. 


(LALL F££ 


PCF 


{£?e>C!>. 


VA-lL'f 


D P F 


• IOOC7. 


pRlCt or FUEL 


PRICE • 


I^S 


Ak/k/uAl i\isura^ce rate 


AXR 


.1 


ASSET LfFt 


NH 


\ 


SALVA&Cr VALUE 


SV 


O 


■ PI5COWMT FACTOR 


D.F 


,Zo 


SuBSv'bV r^ATf 


SUS IDR 


.Ff 


Pt-^c-evr F(^44/Ve^■b 


PF 


,7F 


lA4Tt-RE5T RATE 


RATE 


.JO 


Coat. Tax rate 


CTR 


-He 


INVEST ME/jr tax CAtbIT 


A 1 . 


..on 








Cost w<Ti&HTt'^(r factors 






F^CTOy^ F(P/? O^TF(7 COST 


CCo 




Fa<.toA. foA steel Cost 


CCS 


1. 


FAOTOA FO/e. MA<U4iA/e/CT c. 




( . 








PROGRAM parameters 






OTJjFcnoc StOic.T<OA/ 


|C?P)7 


1 


Ci>JFr<Ai^T' sf A/S tTfu iry 


DELtA(X) 


Act 



Table 16 (cent.) 
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GEORGE G. SHARP DESIGN 
PROGRAM INPUT PARAMETERS 



PARAMETER. 


.SYMBOL 


VALUE 


NAVIMTIOM CoW5TAAIA/T5 






MAX. LExJ&TH 


LEa/ GXM 


f S' o 


MAX. BEAM 


Beam m 


1 oo 


MAX. DM FT 


DR A FXM 










ECOWOMiC- PA'PAMtnfRS 






fp,fic-ht rate 


FR 


5“o. 


GALL F££ 


?CF 


( OOo. 


PA-iLy 7><?AT F£E 
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USE OF THE SUIviT PROGRAI.'S 



Two programs were developed to obtain the output used 
during the analysis stage. The first utilizes the SUI-IT rou- 
tine to develop attractive alternatives. These alternatives 
are identified by the values of the decision variables and 
the objective function. In addition, a second program is 
available which gives a detailed breakdown of a single de- 
sign alternative. This second program does not employ any 
optimization technique. It merely identifies a design in 
terms of its weights and costs. The inputs for each of 
these programs are outlined in the next section. 

Input ^ormat 

The input deck for the STOiT program contains two 
major elements. The first lists the parameters used by the 
optimization technique and the second contains the data used 
by the design model in defining the objective surface. The 
cards are listed in the order of appearance for the input 
deck of the SUI.iT program. (See Figure 25 ) The second pro- 
gram uses only the design data as its input. See Figure 26, 

CARD Is Problem #1 Parameter Card (See Guide to SUTIT) 
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CARD 2 



CARD 3 



CARD 4 



CARD 5 



CARD 6 



CARD 7 



CARD 8 



CARD 9 



! Initial X_ Vector of Decision Variables 
Format (6E12.6) 

Displacement} Length} Beam} Draft} Depth 
Prismatic Coefficient, 

s Initial Xq Vector (Continued) 

Format (6E12.6) 

Velocity. 

I First Design Card 

Format (5- 10. 2, 5X, 12) 

Propulsive Coefficient} Endurance} Distance 
Between Forts; Maximum Shaft Horsepower of a 
Single Engine } Average Fort Delay Time; Maxi- 
mum Stacking of Containers. 

Second Design Data Card 
Format ( 5F10 . 2, 5X, 12) 

Maximum Payload V/eight; Average Weight of a 
Container; Effective Container Length; Effec- 
tive Container V/idth; Effective Container 
Depth; Maximum Containers Stacked Above Deck. 

Third Design Data Card 
Format (5F10. 2,5X,I2) 

Port Call "^ee; Daily Port Fee; Freight Rate; 
Price of Fuel; Annual Insurance Rate; Indi- 
cator of Objective, 

Fourth Design Data Card 
Format (5F10. 2,5X, 12) 

Percent Financed; Subsidy Rate; Interest Rate 
on Loan; Salvage Value of Ship; Discount Fac- 
tor; Asset Life, 

Fifth Design Data Card 
Format (5F 10. 2,5X, 12) 

Investment Tax Credit. 

Sixth Design Data Card 
Format (5F10.2.5X,I2) 

Maximum Ship Length; Maximum Ship Beam; Maxi- 
mum Ship Draft. 
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CARD 10: Seventh Design Data Card 

Format (5F10. 2, 5X, 12) 

Double Bottom Thickness? Center of Gravity 
of Ballast Weight Above Keel? Ship Avail- 
ability? Annual Corporate Tax Rate. 

CARD 11: Eighth Design Data Card, Sensitivity Data 

Format (8? 8. 2) 

Betal? Beta2? V/WO? WWS? V/WT-l? CCO ? CCS; CCM. 

CARD 12: Ninth Design Data Card 

Format (8F8.2) , 

V/XX; XKC-X; V/PI.IIN. 



CARD 13s Tenth Design Data Card 
Format (313) 

NSHAFT? NENGU? REXTRA. 

CARD Iki Eleventh Design Data Card, Sensitivity Data 
Format (14?5.2) 

(DELTA(1),I = 1,14) 

CARD 15* First Option Card (See Guide to SWIT) 

CARD l6: Tolerance Card (See Guide to SUIuT) 

CARD 1?! Second Option Card (See Guide to SUI<IT) 

CARD 18: Problem #2 Parameter Card for Next Program, 

Repeat Other Cards As Necessary. 
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PRO (SRAM DATA DECK structure 
F igure 24 
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SAKPLE DATA CARDS ?0R SUInT FRCGRAM 
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Sample Output 



Figure 27 shows the output to the program giving the 

detailed design, 

1. Values of the decision variables. 

2. Weight break down of the design, 

3. Container configuration. 

4. Cash flow statement. 

5. Value of objective functions. 

Figure 28 shows the output to the SUMT Program. The 

following annotations apply: 

1. The values of the SUMT parameters used for this prob- 
problem. 

2. The initial Xq vector of decision variables with the 
associated objective value (F), Also an internal clock 
is referenced, (Use of design alternative 4) 

3. The starting feasible point (same as #2 in this problem) 

4. Recording of orthogonal moves made by the routine, 

5. The solution to the first subproblem after 10 moves 
r = 1.0. 

6. The solution to the second subproblem after 12 moves 
r = 1./20 

7. Final solution of the first problem. 
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Figure 28 (cent.) 
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Figure 28 (cent, ) 
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Appendix E 

SOURCE LISTING 0? CCMFUTER PROQR.^J'iS 



l6l 



PROGRAM MAIN 



Routine that calls design model and 
outputs a detailed ship design. 

Calls subroutines; 

READIN 

RESTNT 
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